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Abstract. We show the existence of complete, asymptotically flat Cauchy initial data for 
the vacuum Einstein field equations, free of trapped surfaces, whose future development 
must admit a trapped surface. Moreover, the datum is exactly a constant time slice in 
Minkowski space-time inside and exactly a constant time slice in Kerr space-time outside. 

The proof makes use of the full strength of Christodoulou's work on the dynamical 
formation of black holes and Corvino-Schoen's work on the constructions of initial data 
set. 



1. Introduction 

1.1. Earlier Works. Black holes are the central objects of study in general relativity. 
The presence of a black hole is usually detected through the existence of a trapped sur- 
face, namely a two dimensional space-like sphere whose outgoing and incoming expansions 
are negative. The celebrated Penrose singularity theorem states that under suitable as- 
sumptions, if the space-time has a trapped surface, then the space-time is future causally 
geodesically incomplete, i.e., it must be singular, at least in some weak sense. On the other 
hand, the weak cosmic censorship conjecture (WCC) asserts that singularities have to be 
hidden from an observer at infinity by the event horizon of a black hole. Thus, assuming 
WCC, the theorem of Penrose predicts the existence of black holes, via the exhibition of 
a trapped surface. This is precisely the reason why the trapped surfaces are intimately 
related to the understanding of the mechanism of gravitational collapse. 

In [3], Christodoulou discovered a remarkable mechanism responsible for the dynamical 
formation of trapped surfaces. He proved that a trapped surface can form, even in vac- 
uum space-time, from completely dispersed initial configurations (i.e., free of any trapped 
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surfaces) and purely by the focusing effect of gravitational waves. Christodoulou described 
an open set of initial data (so called short pulse ansatz) on an outgoing null hypersur- 
faces without any symmetry assumptions. Based on the techniques he and Klainerman 
developed in [5] proving the global stability of the Minkowski space-time, he managed to 
understand the cancelations among the different components of connection and curvature. 
This eventually enabled him to obtain a complete picture of how the various geometric 
quantities propagate along the evolution. Christodoulou also proved a version of the above 
result for data prescribed at the past null infinity. He showed that strongly focused gravi- 
tational waves, coming in from past null infinity, lead to a trapped surface. More precisely, 
he showed that if the incoming energy per unit solid angle in each direction in an advanced 
time interval [0, 5} at null infinity is sufficiently large (and yet sufficiently dispersed so that 
no trapped surfaces are present) the focusing effect will lead to gravitational collapses. 
Besides its important physical significance, from the point of view of partial differential 
equations, it establishes the first result on the long time dynamics in general relativity for 
general initial data which are not necessarily close to the Minkowski space-time. 

In [p2], Klainerman and Rodnianski extended Christodoulou's result. They introduced 
a parabolic scaling and studied a broader class of initial data. The new scaling allowed 
them to capture the hidden smallness of the nonlinear interactions in Einstein equations. 
Another key observation in their paper is that, if one enlarges the admissible set of initial 
conditions, the corresponding propagation estimates are much easier to derive. Based on 
this idea, they gave a significant simplification of Christodoulou's result. At the same time, 
their relaxation of the propagation estimates were just enough to guarantee the formation 
of a trapped surface. Based on the geometric sharp trace theorems, which they have intro- 
duced earlier in [li] and applied to local well posedness for vacuum Einstein equations in 
[7], [8], [9] and [10] , they were also able to reduce the number of derivatives needed in the 
estimates from two derivatives on curvature (in Christodoulou's proof) to just one. The 
price for such a simpler proof, with a larger set of data, is that the natural propagation 
estimates, consistent with the new scaling, are weaker than those of Christodoulou's. Nev- 
ertheless, once the main existence results are obtained, improved propagation estimates 
can be derived by assuming more conditions on the initial data, such as those consistent 
with Christodoulou's assumptions. This procedure allows one to recover Christodoulou's 
stronger results in a straightforward manner. We remark that, from a purely analytic 
point of view, the main difficulty of all the aforementioned results on dynamical formation 
of black holes is the proof of the long time existence results. The work [12] overcomes this 
difficulty in an elegant way[] 

When matter fields are present, some black- hole formation results have been established 
much earlier, under additional symmetry assumptions. The most important such work 



We remark that Klainerman and Rodnianski only considered the problem on a finite region. We also 
note that the problem from past null infinity has also been studied by Reiterer and Trubowitz 15 , by a 
different approach. 
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is by Christodoulou (see [3]). He studied the evolutionary formation of singularities for 
the Einstein-scalar field equations, under the assumption of spherical symmetry (notice 
that according to Birkhoff theorem, this assumption is not interesting in vacuum). The 
incoming energy of the scalar field was the main factor leading to the gravitational collapse. 
We note also that the work [3j provides a much more precise picture for the large scale 
structure of the space-time, than that available in the general case. 



In all the aforementioned works, the initial data are prescribed on null hypersurfaces. It 
is however natural to study the question of formation of trapped surfaces for Cauchy initial 
data. We first recall that Cauchy data, i.e., the data defined on a space- like hypersurface, 
must satisfy a system of partial differential equations, namely, the constraint equations 
(1.1). The main advantage of using characteristic initial data is that one has complete 
freedom in specifying data without any constraint. Though Christodoulou's results for ini- 
tial data prescribed at past null infinity predicts, indirectly, the existence of asymptotically 
flat Cauchy data, leading to a future trapped surface, it is natural to provide a constructive 
approach to this problem. In this connection, we mention an interesting piece of work |16| 
by Schoen and Yau. They showed that on a space-like hypersurface, there exists a trapped 
surface when the matter field is condensed in a small region, see also [T7j for an improve- 
ment. Their proof analyzed the constraint equations and made use of their earlier work 
on positive mass theorem (especially the resolution of Jang's equation). We remark that 
their work is not evolutionary and matter fields are essential to the existence of trapped 
surfaces. 



The goal of the current paper is to exhibit Cauchy initial data for vacuum Einstein field 
equations with a precise asymptotic behavior at space-like infinity, free of trapped surfaces, 
which lead to trapped surfaces along the evolution. We give the precise statement of the 
result in next section. 



1.2. Main Result. Let E be a three dimensional differentiable manifold diffeomorphic 
to M 3 and (x±, X2, X3) be the standard coordinate system. We also use |x| to denote the 
usual radius function. Let ro > 1 be a given number, 5 > a small positive number 
and £0 > another small positive number. We divide £ into four concentric regions 
£ = E M U s c* U s s U ^k, where 

Em = {x I \x\ < r }, Ec = {x \ r < \x\ < n}, 
Eg = {x I r\ < \x\ < r2}, and £^ = {x \ \x\ > r2}. 

The numbers r\, r2 will be fixed in the sequel such that r\ — ro = 0(5) and T2 — r\ = O(eo). 

A Cauchy initial datum for vacuum Einstein field equations on E consists of a Riemann- 
ian metric g and a symmetric two tensor k (as the second fundamental form) subject to 
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the following constraint equations: 

R(g)-\k\ 2 + h 2 = 0, 



(1.1) 

where R(g) is the scalar curvature of the metric g and h is the mean curvature. 



divgk — d h = 0, 



In order to state the main theorem, we also need to specify a mass parameter mo > 
which will be defined in an explicit way in the course of the proof. Our main result is as 
follows: 

Main Theorem. For any sufficiently small e > 0, there is a Riemannian metric g and a 



symmetric two tensor k on E satisfying (1.1), such that 



1. Y>m is a constant time slice in Minkowski space-time, in fact, (g,k) = (5ij,0); 

2. T,pc is isometric to a constant time slice all the way up to space-like infinity in a Kerr 
space-time with mass m and angular momentum a. Moreover, \m — mo| + |a| < e; 

3. S is free of trapped surfaces; 

4. Trapped surfaces will form in the future domain of dependence o/E. 

Remark 1. The mass parameter m$ reflects the amount of incoming gravitational energy 
that we inject into the Minkowski space-time through an outgoing null hypersurface. It can 
he computed explicitly from the initial conditions as follows (see Section^for definitions): 

1 f 5 

m ° = 4 J l n o| 2 |x(^M ,6>)|^du. 

Remark 2. The regions Sc and Eg do not appear in the theorem. In the proof, we shall 
see that Eg will be e-close to a constant time slice in a Schwarzschild space-time with mass 
too and Ec will be constructed from Christodoulou 's short pulse ansatz. 
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We now sketch the heuristics of the proof with the help the above picture. We start 
by describing initial data on a truncated null cone (represented by the outgoing segment 
between the bottom vertex of the cone and the point u = 1 + 5 in the picture) . 

The part of the cone between the vertex and u = will be a light cone in Minkowski 
space-time, thus its future development, i.e., the white region in the picture, is flat. 

For data between u = and u = S, we use Christodoulou's initial data (from which 
a trapped surface will appear on the top of the light grey region) with an additional 
condition. This condition requires the incoming energy through this part of the cone 
be spherically symmetric. The main consequence of imposing this condition is that the 
incoming hypersurface from u = 5 (represented by segment between the light grey region 
and the grey region in the picture) will be close to an incoming null cone in Schwarzschild 
space-time. 

For data between u = 5 and u = 1 + 5, we require that its shear be identically zero. 
Together with the data on u = 5, one can show that this part of the data will be close to 
an outgoing null cone in Schwarzschild space-time. 

We construct the grey region (pictured above) by means of solving the vacuum Einstein 
equation using the initial data described in the previous two paragraphs. Thanks to the 
consequences of our additional condition, this region is close to a region in Schwarzschild 
space-time. 

We then can choose a smooth space-like hypersurface such that it coincides with a 
constant time slice in Minkowski space-time in the white region (noted as Ejvf m the 
picture) and it is also close to a constant time slice in Schwarzschild space-time (noted 
as the gluing region in the picture). Due to the closeness to Schwarzschild slice, we can 
do explicit computations to understand the obstruction space needed in Corvino-Schoen 
construction and we can eventually attach a Kerr slice to this region. 

1.3. Comments on the Proof. We would like to address now the motivations for and 
difficulties in the proof. 

As we stated in the Main Theorem, the ultimate goal is to obtain a Kerr slice To 
this end, we would like to use a gluing construction for the constraint equations due to 
Corvino and Schoen [2]. Roughly speaking, if two initial data sets are close to each other, 
this construction allows us to patch one to another along a gluing region without changing 
the data outside. Similar results have also been obtained by Chrusciel and Delay in [6]. 
We also note that in [6] , there is a gluing construction with background metric close to the 
Minkowski metric, while in our work the background will be close to a Schwarzschild slice 
instead. J2] extended an earlier work of Corvino pQ which proved a parallel result for time 
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symmetric data. Both constructions in [JJ and [2] relied on the study of local deformations 
of the constraint maps (see Section 2.3). We will also have to make use of the deformation 
technique. Morally speaking, three ingredients are required to implement this approach: 



A reasonable amount of differentiability, say some C k,a control on the background 
geometry. This is necessary since the local deformation techniques can only be 
proved for a relatively regular class of data, see Section 2.3 or [2] for details. 
Some precise information on the background metric. This is essential since even- 
tually we can only glue two data sets close to each other. A priori, one of them is 
a Kerr slice. Therefore, we expect to contract some space-time close to some Kerr 
slice on a given region. We call this requirement the smallness condition since it 
will be captured eventually by the smallness of e. 

The gluing region must have a fixed size. This is important because the local 
deformation techniques work only on a fixed open set. If the gluing region shrinks 
to zero, it is not clear how one can proceed. 



The first ingredient on the differentiability motivates us to derive the higher order energy 
estimates in Section [3j In Christodoulou's work [3], he obtained energy estimates up to two 
derivatives on the curvature components. He also had energy estimates for higher order 
derivatives, but the bounds were far from sharp. For our purposes, we require not only 
differentiabilities but also better bounds. 

The second ingredient turns out to be the key of the entire proof. First of all, we remark 
that the higher order energy estimates are also important for the second ingredient. To 
clarify this connection, we review a technical part of Christodoulou's proof in [4|. One 
of the most difficult estimates is the control on ||V§?7||l°o. He used two mass aspect 
functions fi and [i coupled with rj and rj in Hodge systems. The procedure allowed him to 
gain one derivative from elliptic estimates to close the bootstrap argument. However, he 
had to incur a loss of smallness on the third derivatives of ry. We remark that the loss only 
happened to the top order derivatives. In fact, for lower order derivatives, we can integrate 
some propagation equations to obtain smallness (but this integration loses one derivative 
instead!). In other words, we can afford one more derivative in exchange for the smallness. 
From the above, it is clear that once we have higher order energy estimates, we can expect 
better controls. 

To obtain the second ingredient, we will show that a carefully designed data will evolve 
to a region that is e-close to some Schwarzschild space-time; it is in this sense that we 
have precise information on the background metric. As a consequence, we will be able to 
write down explicitly a four dimensional space in the kernel of the formal adjoint of the 
linearized constraint maps. This kernel is the obstruction space for the gluing construction, 
and we shall use mass and angular momentum (which is of four dimensions in total) of 
the Kerr family to remove the obstruction. In this connection, we want to point out that 



CAUCHY DATA AND GRAVITATIONAL COLLAPSES 



7 



in the work of Reiterer and Trubowitz [15], they also obtained such a near Schwarzschild 
space-time region (which to the best of our knowledge may not satisfy the requirement in 
our third ingredient). 

To obtain an almost Schwarzschild slab, one needs a careful choice of initial data among 
the short pulse ansatz of Christodoulou. We shall impose the following condition: 



namely, independent of the 6 variable. Heuristically, this says that the total incoming 
gravitational energy is the same for all the spherical directions. We expect that through 
such a pulse of incoming gravitational wave, after an advanced time interval [0,5], we 
can regain certain spherical symmetry. We remark that the above condition is the main 
innovation of the paper and will be crucial to almost all of our estimates. 

Finally, we discuss our third ingredient. One may attempt to carry out the gluing directly 
on the space-time constructed in Christodoulou's work. However, when the parameter 5 
goes to zero, one runs the risk of shrinking the gluing region to zero. To solve this difficulty, 
we will further extend the short pulse ansatz beyond the the advanced time interval [0, 5] 
to [0, 1 + 6]. This extension allows us to further solve Einstein equations to construct a 
space-time slab of a fixed size. 

To close the section, we would like to discuss a way to prove the Main Theorem based 
on the work of Klainerman and Rodnianski [12] instead of the stronger result [1] of 
Christodoulou. Since we have to rely on higher order derivative estimates, we can ac- 
tually start with the existence result in [12] which only has control up to first derivatives 
on curvature. We then can assume more on the data, as we mentioned in the introduction, 
to derive more on the solution. As we will do in sequel, we can do the same induction 
argument on the number of derivatives on curvatures. This will yield the same higher 
order energy estimates and the rest of the proof remains the same. 



2.1. Preliminaries on Geometry. We follow the geometric setup in [3] and use the 
double null foliations for most of the paper. We use T> to denote the underlying space-time 
and use g to denote the background 3+1 dimensional Lorentzian metric. We also use V to 
denote the Levi-Civita connection of the metric g. 

Recall that we have two optical functions u and u defined on T> such that 




2. Preliminaries 



g(Vu,Vu) = g(Vu,Vu) = 0. 
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The space-time T> is spanned by the level sets of u and u. The functions u and u increase 
towards the future. We use C u to denote the outgoing null hypersurfaces generated by the 
level surfaces of u and use C u to denote the incoming null hypersurfaces generated by the 
level surfaces of u. We also use Su,u = C_ u n C u to denote the space-like two sphere. The 

notation c}u ' refers to the part of the incoming cone C u where v! < u < u" and the 

notation c]f ' refers to the part of the outgoing cone C u where u' < u < yf' . 

Following Christodoulou [3], for the initial null hypersurface C uo where uq < —2 is a 
fixed constant, we require that ci" ' ^ is a flat light cone in Minkowski space-time and So, Uo 
is the standard sphere with radius \uq\. Thus, by solving Einstein vacuum equations, we 

know the future domain of dependence of ci" ' ^ is flat. It is equivalent to saying that 
the past of Cq can be isometrically embedded into Minkowski space-time. In particular, 
the incoming null hypersurface Cq coincides with an incoming light cone of the Minkowski 
space-time. In the sequel, on the initial hypersurface C Uo , we shall only specify initial data 

one 11 

We are ready to define various geometric quantities. The positive function Q is defined 
by the formula Q~ 2 = —2g{Wu,Vu). We then define the normalized null pair (e$,e±) 
by es = — 2ilVu and = — 2£lVu. We also need two more null vector fields L = Q,e^ 
and L = Qe^. We remark that the flows generated by L and L preserve the double null 
foliation. On a given two sphere S u u we choose a local orthonormal frame (e\, e-i). We call 
(ei, e2, e3, e&) a null frame. As a convention, throughout the paper, we use capital Latin 
letters A, B,C, ■ ■ ■ to denote an index from 1 to 2, e.g. denotes either e\ or e2] we 
use little Latin letters i,j, k, ■ ■ ■ to denote an index from 1 to 3. Repeated indices should 
always be understood as summations. 

Let (j> to be a tangential tensorfield on D. By definition, <p being tangential means that <fi 
is a priori a tensorfield defined on the space-time T> and all the possible contractions of <fi 
with either or are zeros. We use D(j) and D_4> to denote the projection to S^ u of usual 
Lie derivatives Cl4> an d £l4>- The space-time metric g induces a Riemannian metric ^ on 
S u ,u- We use 9I and y 7 to denote the exterior differential and covariant derivative (with 
respect to ^) on Su, u - 

Let (9 a )a=i,2 be a local coordinate system on the two sphere Sq jU0 . We can extend 
9 A, s to the whole V by first setting L(9 A ) = on C uo , and then setting L{6 ) = on 
V>. Therefore, we obtain a coordinate system (u, u, 6 A ) on T>. In such a coordinate, the 
Lorentzian metric g takes the following form 

g = -2n 2 (du®du + du® du) + ^ AB (dO A - b A du) ® {d6 B - b B du). 

The null vectors L and L can be computed as L = d u and L = d u +b A d e A. By construction, 
we have b (u,uq,0) = 0. In addition, we also require H(u,uq,8) = £1(0, u, 9) = 1. 
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We recall the definitions of null connection coefficients. Roughly speaking, the following 
quantities are Christoffel symbols of V according to the null frame (ei, e 2 , e 3 , e 4 ): 

XAB = g{^AeA,e B ),VA = -^g(V 3 e A ,e 4 ),uJ = ^5<V 4 e 3 ,e 4 ), 

X AB = g(^Ae3,e B ),V A = -^(V 4 e A ,e 3 ),w = ^3(V 3 e 4 ,e 3 ). 

They are all tangential tensorfields. We define x' = ^ _1 X> x' = ^~ 1 X and C = \{v ~ V)- 
The trace of x and X will play an important role in Einstein field equations and they are 

defined by tr% = jj AB XAB an d tr% = ^ AB X AB - We remark that the trace is taken with 
respect to the metric ^ and the indices are raised by ^. By definition, we can check directly 
the following identities filogft = ^(77 + 77), D log £1 = co, Dlogfl = to and Db = 4ri 2 £". 
Here, is the vector field dual to the 1-form (. In the sequel, we will suppress the sign ft 
and use the metric ^ to identify £ and C • 

We can also decompose the curvature tensor into null curvature components: 

a AB = R{eA, e 4 , e B , e 4 ), /3 A 

&ab = R(e A ,e 3 ,e B ,e 3 ),f3_ A 

where / is the volume form on Su,u- 

In order to express the Einstein vacuum equations with respect to a null frame, we have 
to introduce some operators. For a symmetric tangential 2-tensorfield 9, we use 9 and 
ti9 to denote the trace-free part and trace of 9 (with respect to ^). If 9 is trace- free, D9 
and D9 refer to the trace-free part of D9 and D9. Let £ be a tangential 1-form. We 
define some products and operators for later use. For the products, we define (6*1,02) = 
$ AC : $ BL \9i) ab (9 2 )cd and (£1,6) = ^ AB (£i)a(6)b- This also leads to the following 
norms \9\ 2 = (9,9) and |£| 2 = (£,£)■ We then define the contractions (9 ■ £) A = 9 A B £ B , 
(9! • 9 2 ) AB = (9 1 ) A C (9 2 ) CB , 9 X A 9 2 = ^ AC ^ BD (9 1 ) ab (9 2 )cd and = & ® & + 6 ® 

£1 — (£i> £2)^- The Hodge dual for £ is defined by *£ = £ A C £c- For the operators, we define 
dyv£ = f A £ A , cu,r<l£ = £ AB f A £ B and (di^6>) A = f B 9 AB . We finally define a traceless 
operator (f®£) AB = (f^AB + {ft) B A - tyr£t A B- 

For the sake of simplicity, we will use shorthands V and R to denote an arbitrary connec- 
tion coefficient and an arbitrary null curvature component. We also introduce an schematic 
way to write products. Let <\> and tp be arbitrary tangential tensorfields, we also use 4> ■ V* 
to denote an arbitrary contraction of 4> and ip by ^ and This schematic notation only 
captures the quadratic nature of the product and it will be good enough for most of the 
cases when we derive estimates. As an example, the notation T R means a sum of products 
between a connection coefficient and a curvature component. 



= ^R(eA,e4,e 3 ,e 4 ),p = -R(e 3 , e 4 , e 3 , e 4 ), 
= R(e A , e 3 , e 3 , e 4 ), a = ^R(e 3 , e 4 , e^, e B )£ 



AB 
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We use a null frame (ei, e2, e^) to decompose the Einstein vacuum equations Ric(g) = 
into components. This leads to the following null structure equations (where K is the 



Gauss curvature of S U)U ): 








Dx' = 


-a, 




(2.1 


Dtr X ' = 


-^ 2 (tr X ') 2 - fi 2 |xT, 




(2.2 


M = 


—a, 




(2.3 


DXrx = 


-^ 2 (trx') 2 -^ 2 |x'| 2 , 




(2.4 


Dr] = 


fi(x- v-P), 




(2.5 


Dr) = 






(2.6 


Du = 


n 2 (2(r ] ,rf)-\ V \ 2 -p), 




(2.7 


Duj = 


2 (2(r / ,r Z )-H 2 -p), 




(2.8 


K = 


--trxtrx+ -(£x) - P, 




(2.9 


dftx = 


L , w 1 , 
2^ tr X - X • V + 2 tr X " 




(2.10 


d^x' = 


I i 
-dtrx' - X ■ V + -trx'r? - 




(2.11 


cu/l?7 = 


1^ ^ 
o- - -X Ax, 




(z.lz 


cu/It? = 


1 

_<7 + A X, 




(2.13 


om) = 


Q 2 (y7®r] + r}(&r] + ^trxx - 


- 7,trxx), 


(2.14 


D(Qtrx) = 


2 (2d^vr Z +2H 2 -(x,x) 


1 

- 2 tr ^ tr ^ + 2 ^)' 


(2.15 


mm = 


9, 2 (f<3r] + r/tgir/ + ^trxx - 


- ^trxx), 


(2.16 


D(Qtr X ) = 


2 (2d^vr ? + 2|r / | 2 -(x,x) 


1 

- 2 tr ^ tr ^ + 2 P>- 


(2.17 



We also use the null frame to decompose second Bianchi identity V[ a Rbc]de = m to 
components. This leads the following null Bianchi equations, 

Da- ^mrxa + 2wa + 0{-yg/3- (4r7 + C)§/3 + 3xp + 3* X cr} = 0, (2.18) 

Da - ^ntr X a + 2ua + tt{f ®fi + (4r? - C)®/3_ + 3xp - 3*x<r} = 0, (2.19) 
3 

Df3 + -ntvxP - tlx- /3 -u/3 - 0{diva + (r? + 2() ■ a} = 0, (2.20) 
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DJ3 + -SltrxP -nx-/3-uij3 + fi{dyva + (rj - 


-2C) -a} 


= 0, 


(2.21) 


DJ3 + 


I 

-ntix/3 -nx-/3 + uj3- fi{dp + *dcr + 3rjp + 3*7/0" - 


+- n-P] 


= o, 


(2.22) 


D§_ + 


i 

-ntixfi - fix • P + oj/3 + VL{$p - *$a + 3r/p - 3*7/0- - 


-2x-P} 


= 0, 


(2.23) 




Dp + -Mrxp ~ fi{d^/3 + (2r? + C, /3) - 




= o, 


(2.24) 




Dp + ^Utvxp + n{dftp + (2?7 - C, £) - 


\ix, «)} 


= 0, 


(2.25) 




L>a + |ntrx<7 + fi{cu/l/3 + (2t/ + £ - 


2^ A a} 


= 0, 


(2.26) 




Off + ^fitrxcr + fi{cu/l/3 + (2r/ - £ *P) + 


2XAa} 


= 0. 


(2.27) 



To conclude this subsection, we recall how one prescribes characteristic data for Einstein 
vacuum equations on two transversally intersecting null hypersurfaces, say Cu ( f +1 ^ U C 
in our current situation. In general, the initial data given on C uo U C should consist of 
the full metric ^ , torsion £ outgoing expansion trx and incoming expansion trx on the 

intersecting sphere Sq jU0 together with the conformal geometry on ci^'* 5 and C . As we 
observed, the incoming surface C is a fixed cone in the Minkowski space-time, thus ^ u , 
CMX an d trx are already fixed on So,u - Therefore, to specify initial data, we only need 

to specify the conformal geometry on Cuq S+1 ^- We will see how the short pulse data of 
Christodoulou is prescribed in next subsection. 



2.2. The Work of Christodoulou. We first discuss Christodoulou's short pulse ansatz 
presented in [3]. As we mentioned earlier, we need to specify the conformal geometry 

on Ci°/ +1] . Let be the one parameter group generated by L. We can rewrite the 

induced metric $\ S&UQ uniquely as ^|s„,„ = (0|s„,„ o ) 2 ?k,„ , where is a positive 

function, such that the metric 3>^|s Ui „ on So,u has the same volume form as $\s , uo - I n 

this language, we only need to specify ^* u i\s u „ freely on Sq <Uq since it yields the conformal 

geometry on cl° ' 5+1 '. 

Let {(U\, {Of )), (U2, (0% ))}a=1,2 be the two stereographic charts on 5o,« . Thus, the 

I |2 

round metric ^|s , uo is expressed as ($\s o , uo )ab(0) = (1 ^Tjg| 2)2 &AB with 9 = 9\ or 2 and 
\6\ 2 = l^ 1 ) 2 + \9 2 \ 2 . Since we require that $(u) = <&* u </)\s uu on So,u has the same volume 
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form as $\s , uo > that is > 

detQ(u) AB )(e)=det(mAB)(e) 



(l + il^l 2 )^ 
(u) is given by 



Nl 2 , nS \u \ 2 



where m^s takes value in the set of positive definite symmetric matrices with determinant 
one and ipAB takes value in the set of symmetric trace-free matrices. After this reduction, 
to prescribe initial data, we only need to specify a function 

-0 : [0, 6 + 1] x S 0)Uo — > S 2 , (u, 6) !->• exp tPab(u, 0), 

where S2 denotes the set of 2 x 2 symmetric trace- free matrices. 

In Christodoulou's work [4J, he only provided data on Cuf \ i.e. u E [0,(5]. First he 
chose a smooth compactly supported ^-valued function ^0 £ C%°((0, 1) x 5o iUo ). Then he 
called the following specific data 

Pol 

the s/iori pu/se ansatz and he called V'o the seed data. 

Remark 3. A key ingredient for the current work is to give further restriction on the seed 



dataipo (see (2.30) or (2.31)). We shall further extend ifi to the whole region [0,5+1] x <So,uo 
by zero. 

To state the main theorems in [3], especially the energy estimates, we also need to define 
some norms. Let k 6 Z>o be a non-negative integer, on each outgoing cone Cu' S \ we define 

K k (u) = r5|ii|-Hl|<55|«|(HV) fc a||+||53| u | 2 (| u |V) fc ^||+|||u| 3 (HV)*(p,ff) +||<T 1 M 4 (MV) fc £||), 
and 

2 

\u\ 



O k+1 (u) = lup 1 | U | (|«| V) fc+1 x|| + ||h| 2 (|n|V) fe+1 (tr X - — )|| + \u\ 2 (\u\V) k+1 x\ 



+ ||rVl 3 (hlv) fc+1 (trx + ri- r4)ll + IN 2 (Mv) fc+1 fo,*/)ll 



lul |it| 2 ' 



+ ||5-V| 3 (|u|V) fc+1 w|| + ||<52|n|(|u|V) fc+1 w| 
.11 

we define 



where all the norms II • II are taken with respect to II • II r2 ,M°,8h', on each incoming cone CL, 

L (C u ) — 



KkhKu) = \\\u\ U Mu\ 9 2(\u\V) k a\\ L 2 { c u )- 
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Recall that we use r or R to denote an arbitrary connection coefficient or curvature 
component. We now use 1Z k [R](u) , O k+ i{T](u) to represent the corresponding norms of 
this given component. To be more precise, let (ft = T or R, we have 

K k [(ft](u) = 5-"2\ u \^\\5-^M- s ^K\n\f) k ^c u ), 
K k [cft)(u) = \\\u\-h-^M- s ^\W\f) k HLHcj, 

o k+ M(u) = r^|«|- 1 ||r'-WH-W(|«|y) fc+ VllL»(c I .)- 

where the r(cft) or s((ft) can be easily retrieved from the precise definitions for each specific 
component. Similarly, we introduce 

ni_M{u,u) = |u|-l||r^)|«r s W(|«|y) fe -VlU 4( 5 H , u ), 
Kf__ 2 [(ft]{u,u) = \\5-^M- sW (W\f) k - 2 H L -(s^y 

Finally, we define total norms 

K k = supK k (u), O k+1 = supO fc+ i(n), K< k = ^Uj, 0< k = ^Oj. 

j<k j<k 



We state the main result in [3]. Roughly speaking, it asserts that up to two derivatives 
on curvature, all the norms (k < 2) defined above can propagate along the evolution of 
Einstein vacuum equations. 

Christodoulou's Main Estimates (Theorem 16.1 in [3]). If 5 > is sufficiently small 
(depending on the C 8 bound of the seed data ipo), there exists a unique solution for the 
Einstein field equations on a domain M corresponding to < u< 5 and uq < u < — 1 — 5. 
Moreover, the following total norms 

n<2, E< 2 [«], <i, k^, o< 3 , o% 2 , oz, 

are bounded by a constant depending on the C 8 bound of the seed data i/jq. 



The main consequence of this theorem is the dynamical formation of trapped surfaces. 
In [1], Christodoulou showed that if 

uniformly in 6, then a trapped surface forms in the future of C Uo . 



2.3. The Work of Corvino-Schoen. We follow closely the notations used by Corvino 
and Schoen in [2] unless there are conflicts with the current work. We first recall some 
definitions on a given three dimensional space-like slice S in a vacuum space-time. The 
vacuum initial data on £ consist of a Riemannian metric g and a symmetric 2-tensor k 
subject to the following constraints (to distinguish from the notations in four dimensions, 
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we shall always use barred notations in three dimensions, e.g. V denotes the Levi-Civita 
connection associated to g): 

R(g)-\k\ 2 + h 2 = Q, 

divgk — d h = 0, 

where h = tigk = gijk 1 ^ denotes the mean curvature, divg/cj = Vjk^ , and all quantities are 
computed with respect to g. One then rewrites the constraint equations by introducing 
the momentum tensor 7fjj = kij — h ■ g~ij. Let T~L and <3? be the following maps: 

n(g,7t) = R s + ^(tr7t) 2 -\tt\ 2 , 

$(5,7f) = (7£(^,7f),div g 7f). 
The constraints then take the form <&(<?, 7f) = 0. 

We use A4 k ' a (T l ), S k ' a (T,) and X k ' a (Yi) to denote the set of Riemannian metric, sym- 
metric two tensors and vector fields on £ with C k,a regularity respectively. Thus, we 
have 

$ : M k+2 ' a (Z) x S k+2 ' a (Z) -> C k ' a (T) x X k+1 ' a (H). 
The formal L 2 -adjoint operator D<§>*,- -■, of the linearization DQ^ ^ is then given by 

DU\ m {f) = ((Ltf)^ + (tr 3 f • fry - 2^^)/, (tr § 7f • ^ - 2%)/), 

Ddiv^j(X) = ^(£xmj + VkX^ - (XiiVk^j + XjVk^i) (2.29) 

- (V m X k ir km + X k V m 7r mk )gij, -£-x9ij), 

where L^f = — Agf + V 2 / — / • Ric(g) is the formal L 2 -adjoint of the linearization of the 
scalar curvature operator. 

Let Q C S be a given bounded domain with smooth boundary and p be a smooth 
positive function on Q which near £1 decays as a power of distance to the boundary, i.e. 
p ~ d N where d is the distance to the boundary and N will be fixed later (see the theorem 
below). The weighted Holder space C k fl(Q) is defined by the norm H/H^a = \\fp~^ \\c k > a 

in the obvious way; similarly, we can define those spaces for tensors. We are ready to state 
the local deformation theorem which will play a key role in our gluing construction. 

Local Deformation Theorem (Theorem 2 in |2j). Let ( G Cq°(Q) be a bump function 
and (go, no) G M k+ ^ a (T) x 5 fc+3 ' a (S). Then for N sufficiently large, there is an e > 
such that for all pair (v, W) G C k > a (Q) x X k+1 ' a (Cl) with the support of (v, W) - <&(go, n ) 
contained in & and with \\{y,W) — ^(go,no)\\ c k, a ^ x x k+1 - c '(n) ^ e ' there is a pair (g,n) G 

M k+2,a^ x s k + 2 > a (E), such that $>(g,n) - (v,W) G C ■ KerDQ* . in Q and (g,Tt) = 
(go, no) outside Q. Moreover, (g,ir) G Ai k+2,a (T,) x 5 fc+2,Q (S) depends continuously on 
(v,W) - $(go,Tro) G C^l"(J2) x X^' a ((Y). 
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In [2], Corvino and Schoen used this theorem to approximate asymptotically flat initial 
data for the vacuum Einstein field equations by solutions which agree with the original 
data inside a given domain and are identical to that of a suitable Kerr slice. We will use 
this theorem in a similar manner in a different situation to prove our main theorem. 



2.4. The Structure of the Proof. This section is devoted to an outline of the proof. It 
consists of three steps. 



• Step 1 Higher Order Energy Estimates. 



As we mentioned in the introduction, since we would like to use Corvino-Schoen construc- 



tion (see the Local Deformation Theorem in Section 2.3) which requires certain regularity 
of the space-time, we are obliged to derive higher order energy estimates. This is accom- 
plished in Section [3} 

We will derive the estimates on M where < u < 5. Christodoulou's work [3] already 
proved such estimates on the level of two derivatives on the curvature (see the Christodoulou 



Main Estimates in Section 2.2). His estimates are already good enough to construct the 
space-time. Thus, in our case, we can use an induction (on the number of derivatives) 
argument instead a bootstrap argument. 



Although this part should be regarded as the routine proof of the persistence of regularity 
for vacuum Einstein field equations, we would like to emphasize that by affording more 
derivatives, we can gain smallness in 5 compared to Christodoulou's work. For example, 
in in L°° norm, we have 

\ u \ 

while, in the current work, since we can control the third derivatives on the curvature 
components, we actually have 

\fv\ < ri3^- 
\ u \ 

The gain in 5 will play a crucial role in Step 2. 



• Step 2 Construction of a Transition Region Close to the Schwarzschild Space-time 
with Mass mo- 



This step is the main innovation of the paper and it is completed in Section |4j Roughly 
speaking, this step build a bridge from Christodoulou's work to the Corvino-Schoen con- 
struction. By solving the vacuum Einstein field equations, we can construct a region with 
a fixed size (independent of the small parameter 5) which is close (measured by 5) to a 
region in the domain of outer communication in the Schwarzschild space-time with mass 
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itlq. It is precisely on this region, or more precisely a space- like hypersurface inside this 
region, where we can use the Corvino-Schoen construction. 

As we mentioned in Remark [3] in addition to Christodoulou's short pulse ansatz de- 
scribed in Section 2.2 we need to impose one more condition on the initial data defined 



on Cuq j that is, for all 9, 



ds 



(8,6) 



ds = 16m,Q. 



(2.30) 



This condition is slightly different from the condition proposed in Remark [TJ namely, 



\ u o\ 



\x(u, uo, 



',du = 4m 



o- 



(2.31) 



The advantage of using the former over the latter is that as 5 changes, (2.30) is always 



valid, but (2.31) is not. However, under the condition (2.30), (2.31) is valid up to an error 



of size 5 to some positive power. The condition (2.31) is more physical and easier to use 



In this paper, although we impose data via the seed function i/jq, we shall stick to the the 
latter condition. 



The condition (2.31) is the key ingredient to Step 2 and most of the estimates are 



directly tied to it. The condition (2.31) also has a clear physical interpretation: it requires 



the incoming gravitational energy per solid angle in the advanced time interval [0, 8] to be 
the same for all angles 9. Roughly speaking, we impose certain spherical symmetry on the 
initial data. 



In Section 4.1 , we show that the sphere Ss jUo is 5-close to a given sphere in the Schwarzschild 
space-time with mass mo; here, the closeness is measured in C k norms. The proof is based 



on condition (2.31) 



In Section 4.2, we show that the incoming cone C_ s is 5-close to a given incoming cone 



in the Schwarzschild space-time with mass mo- This cone is rooted on the given sphere 



Section 4.1. The proof is based on condition (2.31) as well as the higher order energy 



estimates derived in Section [3j 



In Section 



4.3 



we further extend our data on C, 



[8,6+1] 



smoothly by setting 



X = 0, 6<u<6+l. (2.32) 

We then show that the outgoing cone cj^' <5+1 ^ is <5-close to a given outgoing cone in the 
Schwarzschild space-time with mass mo- This cone is also rooted on the given sphere 



Section 4.1. The proof is based on condition (2.31). 



In Section 4.4, since the Sections 4.1, 4.2 and 4.3 provide a characteristic data set for 
the vacuum Einstein field equations, we will solve the field equations to further extend 
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Christodoulou's solution. By virtue of the closeness to the Schwarzschild data, we show 
that the resulting space-time is close to a region in the domain of outer communication in 
the Schwarzschild space-time with mass tuq. Furthermore, we will also see that this region 
has a fixed size, namely, it will not shrink to zero when we decrease the parameter 5. 

• Step 3 Gluing a Kerr Slice. 

This step is completed in Section [5] As a consequence of Section |4.4[ we can choose a 
three dimensional space-like region for gluing. Then we follow the procedure in [I] and 
[2]. Thanks to the closeness derived in Step 2, we can almost explicitly write down the 
kernel of the adjoint of the linearized operator of the constraint map. This allows one to 
use the Local Deformation Theorem in Section [2. 31 The kernel will have four dimensions. 
Combined with a fixed point argument, we use a four parameter family of Kerr space-time, 
namely the mass m and the angular momentum a to kill the kernel. We also show that 
the resulting initial data is free of trapped surfaces. 

3. Higher Order Energy Estimates 

In this section, we derive energy estimates for higher order derivatives on curvature 
components. At the same time, this also yields the control of higher order derivatives on 
connection coefficients. More precisely, we shall prove 

Theorem 1. Let k G Z>3. If 5 > is sufficiently small depending C k+N norm of the seed 
data ipQ for sufficient large N , then the following quantities 

fc<k,fc<k[a],fc<k-i,K< k _ 2 ,0<k+i, 0< fc , 0< fc _i, 

are bounded by a constant depending on C k+N bounds of the seed data. 

The integer N is chosen such that for all k, lZk(uo) is bounded by a constant depending 
on the C k+N bound of the seed data. We remark that for k < 2, the theorem was proved in 
[3] with A = 6. In fact the precise value of A is not important in the current work. Then 
we shall proceed by induction on k. We then make the induction assumption as follows, 
for k > 3 

ft<*-l + K<k-lh) + ^<k-2 + ^<k-3 + 0<k + Ok k -i + 0<fc-2 < iVl+7V, (3.1) 

where Fi is depending only on the C l bound of seed data. 

3.1. Sobolev Inequalities and Elliptic Estimates for Hodge Systems. In this sec- 
tion, we recall the elliptic estimates for Hodge systems on S u ,u- Together with Sobolev 
inequalities, this will serve as the basic tools for us to control the L°° norms. We the refer 
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to [1] for the proof and we shall take 5 > to be sufficiently small as in [4j. For the sake 
of simplicity, we use S to denote Su, u , C to denote C u and C to denote C_ u . 

We first collect all the Sobolev inequalities as follows: 

u\\lhs) s mH?z\\vw + wr^uw^s), (3.2) 
neiu«(s) < nhifchHs) + h^uWlhs), (3-3) 

sup(|^||£|U4 (s) ) < U\\ LHS) + ||^III 2(c) (ll^lli (c) + H 2 \\fC\\l Hc) ), (3.4) 

sM\u\<*u\\l*( S )) < \uo\ q UhHs % , uo ) + \\Hm\\h i0 (\M q ~^\\i i0 + \\W\ q fC\\h { c))- 

(3.5) 

For Hodge systems, we have the following two cases. First of all, if 9 is a traceless 
symmetric two tensor, such that 

dft9 = f, 

where / is a one form on S, thus, 

\\? k+1 o\\ms) < \\f h f\\ms) + lumr'-'fWms) + wmf^w^ (3.6) 

fc-i 

+ \u\- 2 \\f k - l e\\ L , {s) + Y.iWTKh^s) + H-'Wf^KW^Wf"- 1 -^; 

i=i 

Secondly, for 1-form £ satisfies 

= /, cu/l£ = 9, 

we have 

\\y k+1 H\\ms) <\\f k f\\ms) + H-'Wf^fWms) + ll^^ll^(^) + WrHf^gh^s) 

+ KrHiiy'cii^cs) + i^r'iiy^iu^s)) (3.7) 
fc-i 

1=1 

Finally, we also collect two Gronwall type estimates for later uses. They can be found 
in Chapter 4 of @J. 

If we have 

D9 = V -ntr X 9 + 7 • + £, 
with I7I < mfi|x| where ^ and m are constants, then 

ru 

P\\lp{S u ,u) ~W,m II^IIlp(S ,«) + / U\\lp(S u , u )- 

JO 
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If we have 

DO = ^ntrx9 + 7 • + f 

where 6 is a r-covariant tangential tensor field and I7I < mfi|x| with some constants v and 
m, then 

2 2 /* 2 

\ u \ r ~ V ~ 3 W\\Lv{Su,u) %P,r,v,m \-Uo\ r ~ V ~* Mlp(Su,u ) + / \v?\ r ~ V ~* U\\lp(S u<u ,)- 



3.2. Estimates for Connection Coefficients. 



3.2.1. Estimates for f k+1 x' and f k+1 tTx'. We first consider f k+1 x' and f k+1 trx' and 
we will give all the details for these two terms. In the following sections, since we shall 
proceed in a similar way, we will be sketchy and only emphasize the key points. 



In view of the terms on the right hand side of (2.10), by induction hypothesis (3.1), we 
compute 



\f h (x - v)\\l2(c u ) ^ \\f k ^h^c u )hh- + / ^\W 1 x\\l^s^j\W 

Jo i=i 



k-l-i 



+ \\x'\\L°°\\f k v\\L2(C u ) i$ 8 2 I 



We remark that the factor 8s\u\ comes from the norms || • ||x,2(c u ). Similarly, we have 

\\f k Mv)\\L^)<H- 3 - k -s^\u\. 

If one replaces k by k — 1, it is obvious that the above estimates also hold. We denote 
i = —x' ■ r\ + ^tvx'i] — i.e. those terms in ( |2.10 ), we have 



fi\\i^(c u ) + H~ 1 \\f k l i\\v{c u ) < \\f K P\\L*{c u ) + 5~Hu\~' 6 ~ k 



-3-k 



(3.8) 



In view of ( |3.6[ ) and (3.7), we turn to the estimates on Gauss curvature K. For % 
0,--- ,/c-l, by (pi), 



17**11^(0.) % \Tp\\l^) + VP\xM\l\C u ) + \\T(^X^X)\\l\c^ 



< 



-2-i 



■ 5l Id. 



Therefore, 
fc-i 



2 \u\ 



-2-k 



62 Id. 



i=l 



According to (2.10), (3.6) and the inequality J^\\ ■ [^(5 , ) < £2 [| • ||i2(c„)) we obtain 



x\\ms u ,.j< I \\? k+1 trx'\\L>( S ,) + 6ny' 



\ L 2 [Cu) +5 3| u | 2 fe .%|. 

(3.9) 
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We turn to ( |2.2[ ), by commuting derivatives, we have 

Df k+1 tr X ' =[D,f k ]itT X ' - ntrxf k+1 tvx' - 2Qx®? k+1 % (3.10) 
fc-i fc-i 

- Y,T V(^trx) ^ - ® f^fx' 



8=0 



i=0 



^ V7' ; ddogft® (? i trx0y' l trx + 2y7 J x(»y ?i x) 



i,j,h>0 



In fact, the commutator [D, can be written as 



2-1 



J=0 



so 



fc-1 



i=0 



Recall that y*ddogf2 = ^'(77 + 77), we rewrite (3.10) as 

Df k+1 tr X ' = -ntT X f k+1 trx'-2nx^f k+1 x+^o.t. 



(3.1). Therefore, 



The notation l.o.t. means the terms have already been estimated by induction hypothesis 

(3.11) 



\f K ^tTx'\\ L 2 {Saiu) <S-2\ 



\f x IU 2 (s tt / jU ) + M • l«l 



Combining (3.9) and (3.11), we obtain 



||y fe+1 tr X 'lb(S H> „) < 
Thus, 



-i 



y fc+1 tr X '|| L2(5ii , j + \u\- l \\f K P\\l*{C u ) + \u 



l-3-fe 



u\ 



' k+1 trx'hHS^)<H~ 3 ~ k (K k m + C)-\u\ 



(3.12) 



Finally, according to (3.8), (2.10) and (3.6), we obtain 



\n k [p} + c)-5^\u\. 



L 2 [Cu) + 5 l\u\ 2 k 



6? \u\ 



<5-hu\- 2 - k < 



(3.13) 
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3.2.2. Smallness on y7 k (r],rj). In order to f k+1 x and y fc+1 trx', we must improve the 
estimates on \\f k (r], ?7)|| £2 ((?„)• We also mentioned this in the introduction. We can actually 
afford one more derivative to gain this smallness in powers of 5. We also remark that at 
this stage, we can only prove that ??)||l 2 (c„) is bounded without any smallness. 

We commute y fc with (2.5) to derive 

fe-l 

Df k r, = n X f k V + Y^ify (^X)f k ~ 1 ~ i ri + Tf (fix) ® f^v) - ? *W), (3-14) 
Therefore, 



i=0 



1 f— 1 

,o^^ 5 H _1 / iiy*2M^j+* 5 iiyViix a (c„)+%r 3 -*-H- (3-i5) 



Similarly, according to (2.6), we have 



k-l 



Df\ = n K f k V + (VxW^V + Tf (%) ® y fe - 1 - i r 7 ) + y fc (^), (3.16) 



i=0 



which implies 

\k\\Tftk 



u\ k \\Tv\\mc u) <Wo\ k \\f k v\\mc uo) + 



,i\k l./i— lnxftk 



u'n\r K ri\\L>(c u ,) (3-17) 



+ / \u\ k \\? k p\\ L i { c u ,)+ 1 Kr%' 



" fe ^- / '- 4 - fc . < ys|u , |. 



We then substitute (3.17) into (3.15) and take supremum on u. As a result, we have 

supiuhifSiu.^^KI'M-'liySll^^j + ^lnl- 1 / supKl^U^IU^,) 

+ <5i|n|- 4 (^[^] + C) • |u| +(55|n|- 2 (^ fe [/3] + 55|u|- 1 C) • 

Thus, 

M*II7 ^lli^.) < Nl^niySllL^) + ^M" 2 ^ + C) • \u\. 



Recall that rj = —r\ on C UQ . In view of (3.15), we have 

Kl fc ||ySlU*(s a , U0 ) < «*|t*>r 2 (w*[0 + ^Kr'c) • Kl- 

Combining all the estimates above, we finally obtain 



\\f K v\\v(s & , u ), 6-^\f\\\ L2{Cu) <di\u\- 2 - k (n k + c) 



u . 
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3.2.3. Estimates for ^7 k+1 x and ft+i trx'. The estimates will be derived in a similar way 
as for )(7 k+1 x' and y fc+1 tr%' '. By (2.11) and the above improved smallness on J7 k r], we 
have 



,k+l. 



W fe+1 x'llL2(C„) i$ 



' fc+1 trx'|k 2(Cu) + \\f k §\\ L 2 {Cu) + 5\\u\- z - k {lZ k + C) ■ 8*\u\. (3.18) 



We then commute y/ k+1 with (2.4) to derive 

Df k+1 tr X ' = -tttY X f k+1 teX ~ ® ? fc+1 X + (trx) 2 ? log ft + l.o.t. 



where l.o.t. denotes the terms appeared in (3.1) (we shall always understand l.o.t. in this 
way). 

Although y fe+1 logft = lf k {i] + rf) is a lower order term (in derivative), this naive 

estimate will not yield the expected smallness of y/ k+1 tix'- To remedy this defect, we 
introduce an auxiliary bootstrap assumption as follows, 

Auxiliary Bootstrap Assumption: ||y k+l log ft|| L 2( Cu ) < 5\u\~ 3 ~ h A k ■ 5^\u\ (3.19) 

for some large constant A k depending on 7Z k . Therefore, 

|fc+2||x77fc+l + 

,l X \\L 2 (Cu) ~ l u OI IIP L1 A, \\L 2 (Cu ) 

'«o 

,/ifc+2 - 1 ' 



/■it 

iTTrfc+lj. /|| << I I fc+2 1| tT7 /n i / I /|fc+2 cJi /|— 2 nvft fe+lx /|| 

\f * r X Hl2(c„) £ Fo| + II y trx ||i,2 (Cuo ) + / \u\ + -52\u\ \\f tr* Hia(o u ,) 



+ / I'.' 
' iiO 



• *3|«'niy ^iu 2( c7 u ,) + %nc^ + A fc + c) • <53|«|. 

According to (2.15), we have |iio I fc+1 1| y 7 A:+1 trx / 1| (o^q ) ^ 5\uo\~ 3 C ■ ^ 5 |^o|> so 

liy^trx'H^^) < 8\u\- 4 - k {n k + A k + C) ■ 5-2\u\. (3.20) 



Thanks to (|3.18|), we obtain 

Il 2 (CV) 



\f k+1 x l 



< 5?\u\- s - k (n k + C + 6*\u\- 1 A h )-6*\u\. 



(3.21) 



3.2.4. Estimates for y k+1 r] and K+L r). The estimates in this section rely on the bootstrap 



• fc+i 



assumption (3.19) and the estimates established in previous sections. In view of (2.12) and 



(2.13), we have two Hodge systems: 



1 1 

d^?7 = -p + -(x,x) - ft, cuj% = -cr+ ^xAx, 



(3.22) 



(3.23) 



where /U, /i are the mass aspect functions (see [4 J for definitions). To justify the names, 
we notice that the integrals of fi or fi on S uu are the Hawking masses of S uu . In reality, 
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we just take the above equations as definitions of fi and \x. They satisfy the following 
propagation equations: 



Dfi 
Dpt 



-ntrxv - ^tx X p - Jfitrxlxl 2 + \^ti X \ri\ 2 + (2^x ■ V ~ ^x X v) (3-24) 
-fltixp - -OtrxM - i^trxlxl 2 + ^trxjvf + <ty (2^X ■ V - ntrxv). (3.25) 



We commute y with (3.24) to derive 

y fc di^(fitrx) -r?-mrxy fc d^vr/ + l.o.t.. 



Z?y fc /i = - tttrxf k H - ^ti X f V + 2y fc d}v(£lx) • + 2^x • ^ fc d^vr/ 



Since we can control y7 fc di^x) an d g fc diy (Otrx) by (|3_12j) and ( |3.13| ), and we can control 
V7 fc dyvr/ and ^d^vr? by (3.7) and (3.22), we have 



JO 



u . 



Similarly, we commute f with (3.25) to derive 



Df k /i 



fitrxyV- -Otrxy*> + 2y fc d^(ilx) ■ V + 2 ^X' f k ^V 
y fc d^(fitrx) -?7-mrxy fc d^vr/ + l.o.t.. 



In view of fl3.7| ), ( |3.20[ ), ( |3.2l[ ) and fl3.23| ), we obtain 

'VIIl 2( c7„) <|nol fc+1 ll?VllL 2 (c uo) + /Vl fe liyVllz^, } 



+ |-u|" 2 (C7 + ^ fc + 52|-u|- 2 A fc ) • 55|n|. 



Since = —2p+ (x, x) — ^ on C Uo , (3.26) yields 

liyVlb(S s ,„ ) < Inol^-^C + ^lnor 1 ^) • kl- 



Together with (3.26) and (3.27), we have 



Once again, we use (3.7), (3.22) and (3.23) to derive the final estimates 

.5 , 



\f + V\\mc u ) <\ur- k (C + K k + 5*\u\- 6 A k ) • 8* \u\, 
\y k+1 v\\L*(c u ) <\ur 3 - k (C + K k + 6l\u\- 2 A k ) ■ S l 2 \u\. 



(3.26) 



(3.27) 



(3.28) 
(3.29) 
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'fc+1, . „„j xtrk+l. 



3.2.5. Estimates for cu and ^7 +1 uj. Those are the last two connection coefficients. 
We will first establish the estimates for y +1 w and close the auxiliary bootstrap argument 



based on (|3.19|) in this section. We first introduce two auxiliary quantities k and k as 

(3.30) 



follows (see [I]) 

fw = K + d^(^), Auj = k - dy(r(Q/3). 
They satisfy the following propagation equations: 



Dk + Otrx« = -2(!)£fw) + m, 
Dk + ntrxK = -2({lx,f 2 uj) + m. 



(3.31) 
(3.32) 



where 



m 



- 2(d^ (nx),iu) + -d# (ntr X ■ W) - (m 2 ),M - (m 2 ), - pfi (^ 2 ) 

+ A (n 2 (2( V , n) - \n\ 2 )) + dft (n 2 ( X ■ P - 2 X ■ + 3 W - 3 V)), 



and 



1 



m = - 2(d^ (Six), fa) ~ ^ (^ tr X • OjS) - W(fi 2 ), $Jp) + (^(O 2 ), *da) - ^ (fi a ) 



+ f (17 2 (2(r/, g) - |r?| 2 )) - d^ (n 2 (x ■ P - 2 X ■ £ - 3 W - 3 V))- 
We commute with ( |3.31 ) to derive 

-2^x-y fc+ Wy fc ~Vi+i.o.t.. 



For y k+1 to, it is controlled by the definition of k and elliptic estimates (3.7). For y 



rjfe+1 



m, 



the highest order terms are ^/ k+1 m are J jP k f3, J7 k f3, J7 k p, J/ k a and J7 k r/, J7 k+1 r]. They 
can be controlled either by (3.1) , (3.28) and ( |3.29| ). Therefore 



f^tkWviS^) < s\u\-^(c + TZ k + 6 2 \u\-*A k ) ■ \u\. 



\-4-k 



By the definition of k and standard elliptic estimates (3.7), we obtain 

- 4 ~ k (C + K k + 5 2 \u\- 3 A k ) ■ d^ul. 



(3.33) 



by commuting with y/ k+1 ; we have 

Df k+1 \ og n = f k+1 u + Y,f l f (%) • f fc_1 ~ V log 



We now close the bootstrap assumption (3.19) using this bound. Recall that to = Dlog Q, 



fc-i 



i=0 



By integrating, we derive 

\u\ k \\f k+1 logO|| L2(Cu) < C'S\u\- 3 (C + U k + 8 2 \u\- 3 A k ) ■ 5?\v\. 
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We then A k = AC'{C + K k ) in (3.19) and choose 5 sufficiently small such that 4C'5 2 < 1, 



then ||y k+1 log n||£,2( Cu ) < \8\u\ 3 k A k ■ 5z\u\. Hence, the bootstrap argument is closed. 
By the choice of A k , we rewrite KM, dOTb, d3l28b, pl29l and (13331 as 



S— ll |3||v7fc+U 'II i r— A| |2||T?7fe+lW|| , I |2 II Tfr fc+1 

6 |u| || V + trx llrarc,,) + <* 2 \ u \ \\f X\\tf(C u ) + m \\f 



+ u 



T fc+1 



I i f— 1| iSllT^fc+l I 



< 



l-Jfe-1 



^IIl 2 ^) 



(C + ft fc )-tf3|li|. 



(3.34) 



Similarly (we now have (3.34) at our disposal), we also obtain the estimates for 



fc+i 



^+b.„ .. . <h w 



|y w IIl 2 (c u ) 



-2-A 



(C + ^ fc ) -<55|u|. 



3.3. Estimates for Deformation tensors. Following [1], we define the rotation vector 
fields Oi (i = 1,2,3) to satisfy DOi = on C UQ and DOi = on M, and Oi|s ,« the 
standard rotation vector field on the round S 2 . For the sake of simplicity, we will suppress 
the lower index i in sequel. The deformation tensor is defined by = Cog- We also 
define its null components ^°Vab = ^^ab, and Z A = Q$ AB ^tt^b- m particular, Z = 
on C no by construction. For tensor field 0, ^ G </> is defined as the projection of to 
By definition, on C U() , 

Dtr(°V = 20(fitrx), D (0) # - fttrx (0) ? = -Oxtr^ty + 2^ (^x)- 
We also have on P: 

Dtr(°V = 2O(0trx), DZ = -4£ (n 2 (), D^j - fltrx} ^ = -ttxt^H + Ho( n x)- 
The purpose of the section is to prove the following estimates : 

\\f k+1 0\\ < \u\~ k (C + 5 1 2\u\- 2 TZ k ) -6*\u\, (3.35) 
8- 1 \u\ 2 \\f k ix^\\ +6-?\u\\\f k (°% + M||V7 fe Z|| < \u\- k {C + K k ) ■ 8*\u\. (3.36) 
where all the norms are take with respect to || • Hxafc^)- 

We first prove by induction that, for i < k — 1, the following quantities 
\\f i+l O\\ L , {C ^8-'\u\ 2 \\^ 

are bounded by |u| _l C • 6? |u| where C = C{F k -\+N)- The case for i = 0, 1,2 have been 
proved in [3]. We assume the case for all j < i— 1 is true. We commute ^ * with the above 
equations to derive (the last three equations only hold on C UQ ) 

Df i+1 0=f i+1 {ttx)-0 + l.o.t., (3.37) 

DfHx^f =2f i £ Q (ntrx) + l.o.t., (3.38) 

Df^i - ntrxf i(0) f = - <g> f *trf°V + 2f i £ (nx) + 1-o.t., (3.39) 

Dy*Z = -4V7^ (i2 2 C) + l.o.t.. (3.40) 
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Df l+1 o =f l+ \n x ) • O + l.O.t., , 
Df'ti^f =2f i j£ (ntT X ) + l.o.t., 



(3.41) 
(3.42) 
(3.43) 



In the current situation, l.o.t. consists of the terms containing at most (i — l)th order 
derivatives of deformation tensors which are bounded by induction hypothesis, and the ith 



order derivatives of connection coefficients which are bounded by (3.1). Hence, by (3.37) 
and Q3.41P , we conclude \\f i+1 0\\ L 2 (Cu) < |u| _i C • 6?\u\. 



For 



b = fitrx, ^trx, Q X , ^ 2 C> we have 
i-1 

f%</> = £ T<P + E f ] f [0) i ■ T' 1 -^ = toV* + T +1 o ■ 4> + i-o.t. 



3=0 



Since \4. f l ^\<\0\\f 



tion hypothesis (3.1) and the estimates for J7 k+1 0, we conclude 



fO\\f % (/>\, combined with (3.38), (3.42), (3.40), the indue 



s-wwr \^i\\ L , {Cu) , \u\\\tz\\ l , { c k) < wrc ■ 8 



2 \U\. 



Combined with (3.39), (3.43), we also conclude 



For the case when i = k, we can proceed in a same way as to the case when i < k — 1. 



The only difference is that the terms 4> = Qti X , ^tr%, Q X , Q X , ^ 2 C> ^ k+l <t> are not bounded 



by (3.1) but (3.12), (3.13), (3.34). Thus, we have completed the proof of (3.35) and (3.36). 



3.4. Estimates for Curvature Components. In this section, based on the induction 



hypothesis (3.1), we derive energy estimates for curvature components. Together with the 



estimates from previous sections, this will complete the induction argument. Usually, the 
derivation of energy estimates is based on the Bel- Robinson tensors, see [5] or [1]. We shall 
proceed directly by integration by parts without Bel- Robinson tensors. This is similar to 



3.4.1. Energy Inequalities. We commute jfi with (2.20), • 
(2T8]) and pTT9| ) . From $LT8§ and $L22$ , we have 



(2.27) and commute fi with 



E k 3 (a) 



Dfy - n x ■ - mftl a = eUp)- 
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where E k (a) and E k (f3) are error terms which will be expressed explicitly later. Therefore, 
D(\tou\ 2 dn^) +D(2\f Q f3\ 2 d^) = (4d^(nlla-f f3)-inl k a(4lo g n,f ^^ 

+ 2E k (a)) + 2{fy, 2E k (P)))d^. 

We can integrate this identity on M to derive 



\£ a\ 2 + / 2|4/3| 2 < / \£ oa \ 2 (3.44) 

C u J C_ Gun 



+ 



J | - 40^a(ddogf}, 4/3) + (?* a, 2E k {a)) + 2(^/3, 2£ 4 fc (/3))|. 



Similarly, from (2.22), (2.24) and (2.26), we have 



2 

which implies 

f \uf\fy\ 2 + f \u\\\f p\ 2 + \f a\ 2 )< I K| 2 |^/3| 2 (3.45) 

+ If |n| 2 (- T 2 T |^/3| 2 + 2^(^log^,4p4/3 + 4 ( T*0) 
JJm \ u \ 

+ (fy,2Elm + (f oP ,2E k (p)) + (f Q a,2E k (a))). 

We remark that \u\ 2 appears as a weight. We deal the remaining equations in a similar 
way. In fact, we have 

Djfcp + + QdfijfeP = El(p), 

Djfea + \ntr X £ k a + Qcufljfep = E k (a), 

Dfop - n x ■ f p + n(if oP - *U k Q a) = E k (P), 

and 

Djfep - Qx ■ foP + ttdftloa = E k (§), 
D$ a - -ntvxioa + SlfQjfep = E k A {a). 

Therefore, we have 

f \ u \\\t k p\ 2 + \f a\ 2 )+ f H 4 |^|</ \u \\\f oP \ 2 + \f o a\ 2 ) (3.46) 
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+ II K| 4 |-^(|4p| 2 + l^| 2 ) + 2fi(^logO,4^-^*^) 

MM \ u \ 

+ (£ k p,2E k 3 (p)) + (J&cr,2El(a)) + (f P,2E k (f3))\- 



and 



I 2\ufW k §\ 2 +( \uf\f a\ 2 <[ 2|no| 6 |^| 2 (3-47) 
+ If \uf\ - ^\f §\ 2 + m o a(ilo g n,f P) + 2(f P,2E*(P)) + (f a,2El(a))\. 

MM \ u \ 



3.4.2. Estimates for Error Terms. We turn to the bound of the error terms appeared in 
the previous sections. We first recall various commutator formulas in [1] and they will be 
used to derive the exact expression of the error terms. For a traceless symmetric two tensor 
9, we have 

1 ^Jif ^FIC 1 

[£ ,d^ = i(°V y A e BC -<ty B { {0) i 9 AC )- 2 ^°W^o) A , 
[lo, D\e = l z e - o) + m {0) i, e\ 
[loMo = - (0) m%e) + m {0) i,e). 

For an one form £, we have 

[£ ,d^]£ = -d^((°)? • - ^tr(°Vd^, [£ ,curl]£ = -itr(°tycu^, 

[i ,r®]t = -2? {0) i ■ i - (0) Mv£ - l( io) f,fm), 

For a function <f>, we have 

W o ,i]<t> = 0, [4, rm = { - 0) i-*i{o i ^), 

it o ,D}0 = z0, [^ O ,m = o- 

We turn to the estimates for E k (a) and E k (f3). The above formulas allow us to write 
fe— l 

E K a ) = ^E^^Prx)?r _, « + t O E 3 (a) + [D,t ]a - (Sl?®fy - fo^f ®P)), 
i=Q 



where £3(0) is from (2.18). For a traceless symmetric two tensor 9, we have 



£ Q 9= Yl f h O-f i2 0---f il O-f il+1 9. 

tiH h 
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Then first term in E^{a) can be written as 

f il O-f i2 0---f ik O-f ik+1 f(ntrx)-f ik+2 a. 

»lH h«fc+2=fc-l 

We bound this term using induction hypothesis. If ik+% = k — 1 or i^. +2 = k — 1, we bound 
the corresponding factors in L 2 (C U ) and the others in \i ij = k — 1 for j = 1, • • • , k, 
ik+i = fc — 2 or ifc + 2 = fc — 2, we bound the corresponding factors in L 4 (S UjU ) (note that 
when k = 3, the case i^+i = ifc+2 = k — 2 can happen) and the others in L°°; for the other 
possible cases, we simply bound all factors in 

For the second term in E k (a), notice that E®(a) = — 2toa— (4r/+C)ig>/3+3x/0+3 *x<7). 
We use r to be an arbitrary connection coefficient and R to be an arbitrary curvature 

0^k 
We can compute j£, E®(a) by the commutator formulas. All the terms can 

be bounded exactly in the same way as for the first term E k {a) except for 

o k - (-2uj_f k a + n(4 V + c) -f k p- srtx-f k p-snx -f k (r). 

This term can only be bounded by IZk- In fact, this collection of terms contains all the 
terms involving kth derivatives of the curvature components while the other terms are in 
the following form 

f h O -f i2 0---f ik O- f ik +^T ■ f ik+2 R. 

UH Hk+2=k,ik+2<k 

which can be easily estimated by induction hypothesis. 

For the third term in E§[a), since [D,/£ ]a = Y^,i=o vM-i tf-o\$'O a )-> according to 

the commutator formulas, this term can be rewritten as 

Y f h o-f i2 o---f ik - 1 o-f ik ^f-f ik+1 (nx)-f ik+2 a. 

iiH hifc+2=fe-l 

The observation is that none of terms contains kth. derivative of the curvature components. 

~ ^k 

Therefore, {D,/jL ]a can also be estimated by induction hypothesis as above. 
For the last term in E%(a), we can write it as 
k k ~ X fc-i-- 

n(f®£ k p-l ¥®p) = nJ2lo \?HoPoP-lo?®M + l -°-t- 

i=0 



Throughout the paper, we shall use T and R as schematic notations. The expression r ■ R denotes a 
sum of products of this form. 
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where l.o.t. denotes the terms with at most fcth derivatives on the curvature components. 
They can be estimated by induction hypothesis. The main term reads as 

f h O ■ f i2 ■ ■ -f^O ■ f ik< -°^ ■ f ik+1 (3. 

iiH Mfc+i=£:,i fc +i fc+ i>l 

Thus, all terms can be estimated by induction hypothesis except for 

k-i . (y7 fc (°)^./3 + (°)^.y7 fc /3) 
which should be bound by lZ k according to the definition of lZ k . 

Putting the estimates together, we obtain 



// 

JJa 



\l k a,2El(a)\ < 5~ l TZ k [a)(C + lZ k ). 

M 



We move on to E%(f3) which can be written as 

k-l 

£ 4 fc (/3) = ^2^0^) ■ fo' 1 '^ + to&M + [D,jfc]P - (Sldiftioa - fe(Od^a)). 

i=0 

Since the estimates can be derived almost in the same way as before, we only emphasize 
the difference. In face, the third term in E\((3) is slightly different from before because 
[D,fi } / 0. This commutator term contains f k Z which should be estimated by lZ k . To 
be more precise, this term can be written as 

i=0 i=0 

^2 f il O---f ik - 1 0-f ih Z-f ik+1 /3. 

iiH Hk+i=k 

All the terms can be bounded by induction hypothesis expect for i k = k and i k+ \ = k. 
Finally, we obtain 



// 

JJM 



\fy, 2E\ (P) | < 5~ 5 K k \J3] (C + ll k + V k+1 [trx] ) • 

We switch to the second group of error terms E k ((3) — j^fioft an d E%(p, a). We have 
E$(P) = ]T Mo^t) ■ to 1 ' 1 ? + toEi(P) 

i=0 

+ [D, f Q ]P - p + *if Q a) - fantfp + »)). 
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It is necessary to observe that we need the term — j^fioP to cancel the term — ^try^/? 

in j^,qE^(I3). The derivation of the estimates is almost the same as before and finally we 
obtain 

\\ M 2 |^,2£f(/3)- <5m k [p](c + n k ). 

JJm \ u \ 

For E 3 (p, a), we need to argue more carefully because it contains some borderline terms 
(meaning that there is not positive power of 5 in the estimates). First of all, we have 

fe-i 

e "(p) = -2 Y,Mo(^x)f ~ 1 ' t p + f E°4( P ) + [D,f ] P + (mftfy - &(CifyP)), 

i=0 

k— 1 

The terms in [D,jfy\p and [Z?,^q]cj do not have the factor f k Z because p and a are 
functions. This is slightly different from previous cases and the estimates are even easier. 
The key different terms are those 4^q{x ' a ) contained in f^E^p) and j^qE2((t). They 
can be bounded by 7Z k [a] and without any gain in 5. Thus, they contribute to borderline 
terms. Finally, we obtain 

// \u\ 2 U k p,2E k 3 (p)) + (f a,2E k 3 (a))\ <K k [p,a](C + K k [a} + sHK k + V k+1 [trx})). 
JJm 

For the remaining error terms, namely E 3 (p,a), E k ((3), E 3 (j3) and E k (a), they can be 
expressed as 

k-l 

E I(P) = ~ 2 E Moi^X^o^P + toEiiP) + [D, fo\P ~ (VdfifoP - &{tify§)), 

i=0 
k-l 

= ~ 2 E Mo(ntrx)t k 1 - l <? + fo E °s(a) + [D, f Q ]a - (Cknfl&P - faOcaflfi)), 

i=0 

+ [d, &\p + (n(if oP - *U k Q a) - f (n(dp - »)), 

k— 1 
k—1 

E k 4 (a) =\Y,Moi^x)lT~ l ^+l k o E l^) + [A?ok+ mHo§_-to(py®p)). 

i=0 
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We can proceed exactly as before to derive 



\unt" oP ,2E$(p) - -f oP ) + (jfcv,2E§(a) - u f Q a)\ < 5 2 K k [p, a] (C + K k + K k [a] ) , 

M 1^1 1^1 



// \u\ 4 \(f P,2E k (P))\ < 6 2 K k [P](C + K k ), 

JJM 

[[ \ufu k p,2El(§) - < 5ln k [§\(c + n k + n k [a]) : 

JJM \ u \ 

[[ \u\ 6 U k a,2E k 4 (a))\ < 5 3 K k [a](C + K k [p,a] + 5^11, + K k la})). 

JJM 



3.5. Completion of the Induction Argument. We first close the induction argument 
on the L 2 level. In fact, based on the estimates derived in the previous sections, we show 
that if 5 is sufficiently small depending on the C k+N bounds of the seed data tpo, then 

K k ,]Z k [a],O k+1 < F k+N 

where F k+ N is a constant depending only on the C k+N bounds of the seed data. 
Recall that for 1-form £, -fo£ = Ei k+1<k f h O---f ik O- f ik+1 £, then 



/ \u\ 2k w k i\ 2 < [ 

J C u J Cu 



+ 1 EM 2 irei 2 - 

Cu i=0 



Similar inequalities hold for traceless symmetric two tensors and functions. 

We suppose that TZ k ,TZ k [a] < G for some large c onsta nt G to b e fixe d later. In v iew o f 
the above estimates, we multiply 5 2 to (Eiib, 5° to (|3.45h, 5- 1 to (|3.46l) and to (|3.47l), 



then substitute the above estimates into ( 3.44 )-( 3.47) and we choose 5 > sufficiently 



small depending on G and F k _i + N (J34J), so that we obtain the following inequalities, 

K k [a} 2 <K k [a\(u ) 2 + 1, 

K k \0] 2 < K k 1(3} (n ) 2 + CK k [p, a] ( 1 + K k [a] ) + 1 , 
K k [p,a} 2 <K k [p,a](u ) 2 + l, 
K k [0\ 2 + ]Z k [a} 2 < TZ k [§\{u ) 2 + CK k [a\(l + K k [p,a)) + 1. 

We substitute the first and the third to the second to derive 

n k [p] 2 < ^ fe [/3](no) 2 + C^K k [p,o]{u^ + 1(1 + V^[a](«o)2 + l) + 1. 
We then substitute the third to the fourth to derive 



KM < C ( l + VK k [ P ,a](u y + l) + Jn k [p](u ) 2 + 1. 
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Thus, we can choose G depending on TZ k (uo) and hence on C k+N bounds of the seed data, 
such that TZ k ,TZ k [a] < G/2. By continuity argument and the estimates for connection 
coefficients in previous sections, we have obtained the estimates for 1Z k , TZ k [a], O k+ \. 

We turn to the following estimates 

7£<fc_i,7£<fc_ 2 >0<fc>C<fc-i < C(F k+N ). 
We first commute f k ~ l with ( |2.20| ), ( |2.24| ), ( |2.26| ) and ( |2.23| ) to derive 

HM- 1 (||^| M |2 ( ^)(| M |v) fc - 1 /3|| i2(Cti) + [IH^^dtilV)*- 1 ^^)^^) 
+ p-V| 4 (5 J D)(|n|V) fc - 1 ^|| L2(Cu) ) < C(F fc+JV ). 



We then apply % Q X f) to ( 2.18 ), to ( 2.22 ), % Q % D_ to ( 2.19 ) and ^ k l D to ( 2.23 ) to 

derive 



*fc-l 



Dtf l DI3 - fix ■ fio~ lD P - fM^ '-Da = £ 
Dfi^Dfi - fix ■ ^ k o X Dfi + Sldft~ k ~' 
'fi^Da - ^ntrxfio + nf® 



tk-l 



-k-l 



1 Df3 


= Et 




l ~ 
Da 


= Et 




l ~ 
Da 


= Et 




1 D(5 


= Et 


- 1)3 («) 



Therefore, we have 



Cy 



\to 



X Da\ 2 < [ \% l Da\ 2 + If \ - AVlf Q 1 Da($ log fl, tfc^DP) (3.48) 

JC un JJM 



-ife-1 . 



+ {f Q ^Da,2Et m (»)) + 2(f -'D^2Er^W))\, 



k-i 



n (fc-l)4. 



and 



u\ s \f - v Da\ z < I 2\uo\*\Po~ l DW + 



fc-i 



o 



16, 



n w ^(^logfi,^- 1 ^) + 2(r - i m 2E r i,o m + {f 'D a ,2E^>°(a))\. 

IM 

This leads to 

5-^\u\\\5^\u\{5D)(\u\f k ' l )a\\ L 2 {Cu) + || j«|-f tf-f (H,D)(|«| V)*- 1 ^!!^^ < C(F fc+iV ). 
The Sobolev inequalities (El, (El and (Ed) yield K^K^ < C{F k+N ). 



M 
k-l 



,(k-l)3. 



(3.49) 

n(fc-l)3. 



Similarly, we commute y ft with ([2JJ), ([2^2], ((^J), ([2J]), (|2_14j) and ( |2.15[ ) to derive 
5-l\u\- 1 (\\5l\u\{8D){\u\y) k x\\^[C u ) + |||n| 2 (^)(|n|V) fc tr X || i2(Cu) 

+ \\6-h\u\ 2 (SD)(\u\V) k x\\ LH c u ) + p-V| 3 (^)(|n|V) fc trxllL 2 (C u ) 
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+ ||M 2 (M))(|u|V) fe 7 / || L2(Cti) + ||rVl 3 (^)(MV)^|| L2(ai) ) < C(F k+N ). 
According to the relation rj + 77 = 2d4og$7, we can also derive 

r*\u\-^\u\\6D)(\u\V) k V \\vr Cu ) < C{F k+N ). 



Now (3.3) and (3.4) yield 0^,0^° < C(F k+ N) except for the component u. 



In order to bound uj, we commute ^J k with (2.8). We also take y fc to both sides of the 
equation and we hope to bound L 2 (C U ) norm for DJ7 u. The commutator [D, y ]u) is easy 
to control thanks to the induction hypothesis (3.1 ). Because we do have L 4 (Su tU ) bound for 



y rj and y rj at the moment, then we can bound || \u | 2 |u| 2 (|u|y ) k ((r], rj) 



< 



C(i*fc + jv)- in view of (3.45), we have |||«| 2 \u\ 2 (\u\y ) k p\\ < C(F k+ j\f) and we can deduce 
that 

\\\u\-l\u\(\u\D)(W\¥) k uj\\ < C(F k+N ), 
Finally, Q and Q yield C^H.O^Jw] < C(F fc+JV ). 

We have completed the proof of Theorem [T] 



4. Construction of the Transition Region 



We have already mentioned that, in addition to Christodoulou's short pulse ansatz 
described, we need to impose one more condition on the initial data defined on cIq'^, that 
is (2.30) or (2.31). This condition is the key to all the construction in sequel. We start 
from the condition (2.30) which is imposed on the seed data, that is, 



Os 



■(8,0) 



ds = 16mo- 



(4.1) 



We show that this condition is equivalent to the more physical one (2.31) up to an error 
of size 5 to some positive power. More precisely, we have 

Lemma 1. On the sphere S^g, f or oil k, we have 

rS 



Kl \x(u,uo)\ du-4m )\ < k 5z\u \ 







The notation A < k B means that A < C k B where C k is a constant depending on the 
(jk+N norms Q f the seed data for sufficient large N. 
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Proof. In stereographic coordinates, we have 



4 / \x(u,uo)\ du = 
10 Jo 



^AC^BD _ 



D$ AB D$ CD du = I m m dy^ABd^CD&u- 
Jo 



In view of the short pulse ansatz (2.28), we have \dg(m AB — 5 AB )\ <fe 5^ \uo\ 1 ■ Now we 
can make use of the definition m = exp ip to derive 

\d$(4 [ S \u \ 2 \x{u,u )\ 2 du- f \u Q \ 2 \d^{u)\ 2 du)\ < k 8*\uo\- x . 
Jo Jo 



In view of (2.28) and (2.30), we have J Q \uQ\ 2 \dui>\ 2 du = \d s ipo\ 2 ds = 16mo, therefore 



\u \ k \f k ( \u \ 2 \x(u,u )\ 2 du-4m )\ < k ^\u \ \ 
Jo 



Thanks to the lemma, we shall not differentiate conditions (2.30) and (2.31). We turn 
to the geometric consequences of these conditions. 



4.1. Geometry on Sg tUQ . The purpose of this section is to prove the following lemma 
which, roughly speaking, says the geometry of the two sphere Sg tUQ is close to the geometry 
of a given two sphere in the Schwarzschild space-time with mass itlq. 



Lemma 2. On the sphere Ss, Uo , we have a = and for all k 



\uq\ 



f (tr*-(i 



\ u o\ 



\ u o\ 



Pol 



\,\uo\\f k o-^J < k <^|u | 



Remark 4. The connection coefficients and the curvature components appearing in the 
lemma have improvement on "smallness" , which means that they are controlled by 5? 
instead of 5 to some nonpositive power in [1] . This improvement comes from the condition 
(2.31). The other connection coefficients and the curvature components, which are not 



mentioned in the lemma, have already been bounded by 5 to some positive power in [3]. 
This suggests that the geometry of the sphere Ss, Uo is close to the geometry of a sphere in 
Schwarzschild spacetime. 



Proof. First of all, a = on Sg Uo follows immediately from (2.1) and the fact that x has 
compact support in C M ° ' S \ 



For tr%, on C Uo , we can write (2.2) as 



Dtix 



\x\ 
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since where £1 = 1. We integrate this equation along [0, 5] to derive 

2 1 f s f 5 

trx(S,u )--, — r = -o/ (tTx{u,uo)fdu- / \x(u,u )\ 2 du. 

Since |ito| fc |y fc (trx— ^ )| |iio| -2 (see Chapter 2 of [4j), taking f k on the above equation 
yields the desired estimates for tr% on the sphere Sg tUo . 

For (3, on C uo , we can write ( |2. 10 ) as 

= ^trx -%•?? + ^trxr/ - /3. 

Since |uo| fc |y fe7 ?| ^Sfc ^ 5 |^o|~ 2 (see Chapter 2 of [4]), x = on S , 5 )Uo , combined with the 
estimates on trx just derived, we can apply f k on the above equation to derived the desired 
estimates for (3 on Ss ;Uo . 



For a, Ton C uo , we can write (2.12) as 

cuflrj = a - Ax- 
Since |uo| fc |y 7A '^| <fc 55|uo|~ 2 and x = 0on <Sv5,u , we obtain the desired estimates for a on 
Sg,u by simply applying f k to the above equation. 



Finally, we consider p. On C Uo , the Gauss equation (2.9) can be written as 

K = - -trxtrx + ^ (x, x) ~ P- 



We rewrite this equation in the following renormalized form 

( K ~ n) + 7 tr x( tr X + r—r) = ~(P ~ Jjj— I fa* - r—r)). 
\Uo\ 4 — | z*o I *\uo\ \uo\ 

Since \u \ k \f k (K - < k S^uo]' 3 and |u | fe |y fe (trx + ^)| < k 5\u \- 2 on C UQ (see 

again Chapter 2 of [4]), combined with the estimates for tr% derived above, we obtain the 
desired estimates for p by applying y k to the above equation. This completes the proof of 
the lemma. I 



4.2. Geometry on C^. We shall study the geometry of C$ by integrating the equations 
( 216) , Q2.17D , Q2.25D , fl2.27| ), ( |2.22[ ), ( |2 . 1 8[ ) , (^ along C 5 . These equations can be viewed 
as ordinary differential equations for the connection coefficients and curvature components 
which are not bounded by 5 to some positive power in [4] . We will see that they are actually 
bounded by 5 to some positive power under the additional condition ( |2.31[ ). We will regain 
the "smallness" from Theorem [I] and improved smallness from Lemma |2| More precisely, 
we prove the following lemma: 
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Lemma 3. Suppose that 5 > has been chosen such that Theorem^ is true for some k. 
Then on the incoming cone C§, we have the following estimates 

\\(\u\f) k - 3 (nx)\\L^<S^\u\-\ fori<k-3, 

l^r l ]|(hiy) fc - 2 (^x)]|^cs 5 ,.) + ^r^Kluiy ) fe - 1 (nx)[k 2( s^) < <5l]nr 2 ; 

\m?nP+¥w,<r)\\L~ <6*\u\- 3 fori<k-3, 
\ u \ 

l«r*ll(My) fr - a (p,^)IU^.o + Itil^lKNy )*- 1 Cp,o-)|| £a(lS r, >0 < ffilt*!- 8 ; 

\\(\u\fyPh°o <6*\u\- 2 fori<k-4, 
\\(\u\jya\\ L oo < 5^\u\- 2 fori<k-5, 

\\(\u\f )*(tr X - (j^i - - < S3\u\~ 3 fori<k-3, 

\u\ \u\ z \u\ z 

\u\-^\\(\u\f) k -\tr X ,co)\\ LHSsu) + |u|- 1 ||(Hy)*- s (tr XlW )|| Ji , (5li<i) <^H- a . 
Remark 5. The smallness gained in the lemma comes from Lemma therefore fi 



rom 



the condition (2.31). Heuristically, the worst control for the geometric quantities on C u 
are propagated from the initial surface C Uo . This is a very important technical feature of 



the work of Christodoulou [3J. Therefore, the condition (2.31) allows one to improve the 
geometry on S$ UQ hence the geometry on C$. 

Proof. We will derive the estimates for Xi Pi °~i ft, a, u and trx one by one. The proof 
for each component is by induction on the number of derivatives. The parameter 5 > is 
chosen to be sufficiently small so that Theorem [T] holds for k 3> 5. 



For x, we rewrite ( 2.16[ ) as 



D(n X ) = ^trx(^x) + 2(Ox, nx)$ + n 2 (f®ri + r/gr? - -ti X x), (4-2) 
where the second term is the trace of D(Q,x)- 

Remark 6. In Christodoulou 's work, we only obtained the bounds for the second derivatives 
of the curvature components which implied the L°° estimates \^ ®rj\ < |u| -3 . If we have 
estimates for the third derivatives of the curvature components, as stated in Theorem 
we have |y g)ry| < 52|^|^ 3 . Therefore, we can bound f2 2 (y ' <8>t/ + 77(8)77 — gtrxx) in L°° by 
5 to some positive power. So the right hand side of the above equation can be treated as a 
small term. This illustrate the idea of the proof and the importance of higher order energy 
estimates. 
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We go back to the proof and integrate D(^x) along C$. The first two terms on the right 
hand side can be absorbed by GronwalPs inequality. In view of the fact that %(<5, uq) = 0, 
according to Theorem [TJ we have 

f u 1 i 

\u\\£lx(S,u)\ < / \u'\\n 2 (f® V + T]®r]--trxx)\du<52\u\- 1 . 

Juo 1 
We then use an induction argument to obtain control for higher order derivatives. Suppose 
that for all j < i < k - 3, |up \\f j (^x)h^(S,u) ^ ^\u\~ 2 , then we apply f i to §t/2§ to 
derive 

+ f i (n 2 (f®r] + r]®r} - ^trxx)) + Lo.t.. 

Therefore, the third term on the right hand side is bounded by \u\ l \^7 1 {Q 2 {^7 ®ri + r]®r] — 
?;trxx))| ^5 <5z |n| -2 thanks to Theorem TJ The l.o.t. is of the form 

j=0 j=0 

We can bound (fitrx) and (fix) by Theoremjljand bound y l ~ J ~ 1 (f2x) by the in- 
duction hypothesis. By integrating the above equation and using the GronwalPs inequality 
to absorb the first two terms, we obtain that 



\u\f)\nx)\\L«<S'\ 



(4.3) 



for i < k — 3. 

For i = k — 2, we only expect that the L 4 (S$ :U ) norm on x is controlled by 5 to some 
positive power, and for i = k — 1, only the L 2 (Ss, u ) norm is small. This is because we 
only have information on IKIwiy 7 )^ -1 ^!!/^^ IKM^ ^^Ki^fS^) < 5? |u|~ 2 from Theorem 
[I] but no information for higher derivatives. By a similar argument and integrating along 
C_$, we obtain 



\u\-^\\(\ u \f) k - 2 (m\LHs s , u ) + M- 1 \\(\«m k - 1 m\\ifip t j z 

It is easy to see that (4.3) and (4.4) also hold if we replace fix by x- 



(4.4) 



We perform a similar argument for the other components. We now bound p and a. We 
recall ^25\) and (pOTl), 



Dp+^ntvxp = -n{dftp + (2 V - c,§) - 2 (x,a)}, 



Da + -fitrxc 



-fl{cufl^ + (277 — (, */3) + ~x A a}. 
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Because k > 3, we have \y j3\ < 5\u\~ 5 . Thus, 

| - n{dftP + (2r/ -(,£)- \{Z Q)}| + I - fi{citfl£ + (27? " C, *£) + ^Aa}|< %|" 5 . 
We renormalize the first equation in the following form 

Dip +^) + %r X (p +^) = ^(Vtrx + n) - + (2t? - C, /?) - J(£ «)}• 

lu^ 2 — | ti | |tj| c> — |it| — — 2 

In view of the facts that |fitrx + A| < <5|7j|~ 2 and the improved smallness for initial data 
on Ss,u in Lemma [2j we integrate the above equations along C_ s and we obtain 

2rriQ 



\P + 



u 



< J 2 1 tj, I 3 and I cr I < <5 2 1 ^ I 



Similarly, we can also run an induction argument as above to obtain 



Hff(p + 



2mo 



< 52 



u 



for i < k — 3, 



It 2 



N?)^ a (p^)||^(fl a ,0,H- 1 [|(Hy)*- 1 (p,a)[| £a(Sji j<tfi|ti 



-3 



(4.5) 



We now bound (3. We examine the terms in (2.22), 



Dfi + ^Otrx/3 - fix • /3 = -ufi + ft{dp + *da + 3t/p + 3 V + 2£ • §}. 



If > 4, we have the bound |dp + *</kr| < #2 |u| 4 thanks to (4.5) and the bound | — u(3 + 



fl(3r]p + 3*7/0" + 2% • /3)| < (52 |?j|~ 5 thanks to Theorem 1 We remark that there is a loss 
of derivatives here because we control (3 in terms of derivatives of p and a via Bianchi 
equations. We then integrate the equation along Cg, combined with the estimates on the 
initial data on Sg )Uo obtained in Lemma [2j we have 



\\(\u\f yf3\\ L ~(S s , u ) < S^\u\- 2 for i < k - 4, 



(4.6) 



Notice that we only have L°° estimates on [k — 4)-th derivatives on (3 instead of (A; — 3)-th 
derivatives due to the derivative loss we just mentioned. 



We now bound a. We rewrite (|2.18|) as follows, 

-n{-f §/3 - (4t/ + 0®0 + 3xp + 3*X°} 



1 



Da — -Otr^a + 2wa — 2(0%, a) 



where 2(fix> Q! ) is the trace of Da. If fc > 5, by (4.6), we have |y<8>/3| < 5^\u\ . We 



remark that at this point there is a loss of derivatives here because we use the derivatives 

. i , 



of (3. We must also use |x| + \(3\ < 52\u\~ 2 by (4.3) and (4.6) to get | - tt{-f®l3 - (4r/ + 
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C)®/? + 3xp + 3*xo-}| < <5 2 — 3 . Therefore, in view of the fact that ce(5, uq) — on Sg iUQ , 
by integrating the equation along Cg, we obtain 

WQulffaU^s^ < ^|n|- 2 for i < k - 5, 

i«r»ii(i«iy ^oiu^j, M^iidtiir)*- 8 "!!^*,.) ^ *»iu|- a . 

We remark that the terms — jfitrxa + 2wa — 2(0%, a)^ on the left hand side are absorbed 
by Gronwall's inequality. We also want to point out that we will not use the L 2 (Sg jU ) 
estimates on y fe ~ 3 a in sequel. Notice also that we only have L°° estimates on {k — 5)-th 
derivatives on a due to the derivative loss we just mentioned. 



Finally, we bound trx and ui. We renormalize (2.17) and (2.8) as follows 



it mr 2 — u it r 



= - ^°)(^tr X + A) + (2f2 2 p + + f} 2 (2d^ + 2|^| 2 - (x, *)), 

\u\ \u\ — \u\ \u\° — 

We have already obtained the following estimates: |fitrx + nr| < £M~ 2 5 |^ 2 /o + xffil ^5 
54- |ii|~ 3 , |n-l| <S\u\- 2 , |0 2 (2d^r ? + 2|r ? | 2 - (x,x)l + l^ 2 (2(??,r/) - |r/| 2 )| < ^\u\~ 3 and 
|fitrx((5, uq) — (|^| — )| < <52|u| -2 . In view of the fact that ito) = 0, by integrating 
the equation along Cg, we obtain 

lltrx-(n-^)lu- + ll--r^ll^^^H" 

\u\ \u\ z \u\ z 

where fitrx has been replaced by trx thanks again to |0 — 1| < 5|ii| -2 . Similarly, we can 
derive estimates for higher order derivatives. Finally, we have 

\\(\u\fY(tr X - (A - ^£),u>- ^)\\ L °o {Ssu) < 5^\u\- 3 for i < k - 3, 

\u\ \u\ z \u\ z ' ' (4 

We complete the proof of the lemma. I 



4.3. Geometry on Cuf +1 \ We extend the initial data on cl° ' 5 ' to cl° ' 5+1 ' by setting 

X = for u E [5, 5 + 1] and we will derive the estimates on Cuf +1 ^ for all the connection 

coefficients and curvature components. For the sake of simplicity, we denote C Uo = Cu<f +1 ^ 
in this subsection. 

We first give another initial data on C Uo U Cg and it is in fact an initial data set for the 
Schwarzschild space-time with mass rriQ. To distinguish this new data from the old one, we 
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shall use a lower index mo for its connection coefficients and curvature components. The 
new initial data set consists of the following quantities: a $ (5,uq) of Sg iUo which is the 
round metric of the sphere with radius \uq\, the torsion Cm (di u o) = on Sg }U0 , two null 
expansions trx mo (<5,uo) = j^y - and trx mo (S,u ) = on S SiUo , the lapse function 

f2 mo = 1 on C Uo U C$ and two shears x m „ =0 on C UQ and x = on C_ 5 . 

We also impose the the condition \uq\ > 2mo- This condition guarantees that the future 
development of the data is in the domain of outer communication of the Schwarzschild 
space-time. We remark that the only nonzero null components for the new data on C UQ 
are trXm , ^ Y X mQ i u mo an d Pm with trx mo > 0, an d the only nonzero null components for 
the new data on C s are trx mo , ^ T X mo > ^m an< ^ Pm ■ 

We now restate the estimates on the sphere Ss, Uo derived Lemma [2] as well as the esti- 
mates given in Chapter 2 of pi]: we have 2 = 0a = 0, w = 0, and for all k, 

iy fc (x,^,^w-w mo ,trx-trx mo ,trx-trx mo ,/3,p-pm ,cr,/3,a)| <k (4.9) 

Notice that we have dropped the weights |«o|- We will not use the weights in sequel since 
we can fix a \uq\ > 2mo now and they are not relevant from now on. We shall prove that 
the above estimates actually hold on C Uo : 

Lemma 4. On the incoming cone Cuf +l \ we have x = a = 0, w = 0, and for all k, we 
have 

\f k {x,il,V,u- io mo , trx - trx mo , trx - try , 0, p - p^ , a, /3, a) | < k S$ . 



Proof. By construction, we also have x = 0, a = and u = on C Uo . The idea to prove 
the estimate is to use null structure equations in a correct order. 



We first control trx- According to (4.9) and choosing 5 sufficiently small, we have 
trx(<5, uo) > so that trx can be solved for u G [(5,(5 + 1] and < trx(:u, uq) < trx(<5, uo). 
We then rewrite (2.2) as 



D (trx - trxmo ) = - - (trx + trx mo ) (trx - trxm ) • 



(4.10) 



Because (trx + trx mo ) is bounded pointwisely, thanks to Gronwall's inequality and ( |4.9[ ), 
we have 

2 4 

|trx-trx mo | < \tTx(S,u ) - (- 



< 55. 



\ u o\ 



\ u o\ 



We can also apply to the above argument and use an induction argument as before to 
derive 



|y fe (trx-trx mo )| < \f k (ti X {S,u ) 



\ u o\ 



A))l<^i 
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We now bound r\ = —rj. We eliminate f3 in (2.5) and (2.10) to derive 



Drj + tixq = -^dtrx- 
Therefore, by Gronwall's inequality and an induction argument, we obtain 

\f k ( V ,rf)(u,uo)\ < k sK 

Before going to the next component, we have to digress to the estimates for the Gauss 
curvature K on Su tUo . Recall that K satisfies the following propagation equations (see 
Chapter 5 of 01): ' 

DK + ntv X K = d^vd^v (Ox) - ^f(Otrx), 



We rewrite the equation on C UQ as 

D(K-K mo )+tr X (K-K mo ) = 
which yields \f k (K(u, u ) - K mo (u,u Q ))\ < k 5^. 



(trx - tr Xm„)K mo - -^trx, 



To bound x and x> we eliminate p in (2.15) and (2.9) to derive the following propagation 
equation for trx — trx 



-m 



D(trx - trx ) + trx(trx - trx ) = 2(K mo - K) - d^?? + |r?| 2 , 



and we rewrite (2.14) as 



These equations yield \ f k (trx(u, u ) - trx m (u,u ))\ + \f k x(u,u )\ < k 6$ . 



We then use (2.10), (2.11) and the estimates above to derive |y (/3, /3)(u, uq)\ <fc 5 2 - 



We rewrite (2.9) as 

K — K 



■-(trxtrx - trXmo tr X mo ) ~(P~ Pm ), 



which implies that \ f (p(u,u ) - p mo (u,u ))\ <k 5 2 - 



We rewrite (2.7) as 



D(u - w ) = -3|r?| 2 - (p - p mo ), 



which implies \ f (u(u, u ) - u m Ju, u ))\ < k Si . 



Finally, (2.19) implies |y a(u, uq)\ < k $ 2 - We complete the proof of the lemma. 
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4.4. Construction of the Transition Region. We consider a characteristic initial data 



problem on Cu<f UC_ 
1- 



data on C^ 



[no-l-6] 



where the data are given by with x = on C, 



[5,5+1] 



_,, . and the 

5 ' is induced by the solution constructed by Christodoulou in [1] under the 
condition (2.30) with |uo| > 2mo. We continue to use the space-time metric g to denote 
the solution of this problem. 



Adapted to this problem, we introduce the following norms: 



= \\T k ( a , 0> P ~ Pm>,<r, P)\\LZ(Ci)> = \\f k (P,p-pm ,<r,P,Qd\\L2(CZ), 



[uo ' u] We take a 
and 



We use short hand notations C« referring to c|f'~' and C" referring to Cu 
small parameter e to be proportional to 5 2 , then the results in Section 



4.1 



Section 



4.2 



Section 4.3 can be summarized in the following proposition: 



Proposition 1. Fix k ^> 5. For any e > 0, there exists 5$ > depending only on the 
fjk+N i) 0un d s of the seed data for some sufficient large N and on e, such that for all 5 < 5q 
and for all (u, u) £ [S, 8 + 1] x [uq, — 1 — S], we have, 

^<ti 2 («o) +^.<lZ S 2 (S) + 0< k - 2 (u,u o ) + 0< k „ 2 (8,u) < e, 

i|TT7<fc— 4 11 , <k 11 ^ 

VP- a \\L*(s S!U ) + \\y- «IIl-(5^ ) < e. 



Based on this proposition, we will prove the following theorem in this section. 

Theorem 2. Let k 3> 5 and e>0. If 5 > is sufficiently small depending on C k+N 
bounds of the seed data for sufficient large N and on e > 0, then there is a unique smooth 
solution (M eo , g) of vacuum Einstein field equations to the characteristic initial data prob- 
lem described above. The space-time M eo corresponds to the region 5 < u < 5 + eo and 
uo < u < uo + £q. The parameter eo > depends only on tuq and uq and it is indepen- 
dent of 5 when 5 is sufficiently small. Moreover, the space-time (M eo ,g) is e-close to the 
Schwarzschild space-time with mass mo in C k ~ 3 norms, namely, 

llff-5m |lc fc - 3 (M 6o , Sroo ) < e - 

Remark 7. We would like to emphasize that once itlq and uq is fixed, the size of eo is 
also fixed. In the introduction, we mentioned that in order to use the Corvino-Schoen 
construction, the gluing region must have a fixed size when 5 goes to zero. The size of eo 
meets this requirement. 

Remark 8. It is not hard to see that g mo is isometric to the Schwarzschild space-time 

with mass niQ. In fact, the initial data for g mo on Cuf +e °^ UCj is spherical symmetric, 
then g mo is also spherical symmetric. The orbits of the isometry group SO (3) are simply 
S u ,u- Therefore by Birkhoff Theorem, g mo is isometric to the Schwarzschild space-time. 
The mass parameter uiq can be figured out by computing the Hawking mass of Ss^ Uo in gmo ■ 
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Proof. The existence of local solutions is due to the classical result [14J of Rendall. We 
make use of a recent improvement [13] due to Luk. Rendall's original proof only gave the 
local existence around the intersection of two initial null hypersurface and Luk showed that 
the existence region can be enlarged to a full neighborhood of initial hypersurfaces in their 
future development. The main advantage of [13] is that it also gives the estimates on the 
solution which depend only on the initial data. We now rephrase the results of [13] in our 
situation: there exists eo which depends only on ijiq and uo and is independent of 5 hence 
e, such that we can solve the metric g for (u, u) £ [5, 5 + eo] x [uo, uq + eo]. Moreover, we 
have the following bounds: 

Proposition 2. We have a smooth solution (M eo ,g) where the M eo corresponds to (u, u) £ 
[5, 5 + 1] x [uq, —1 — 6]. Moreover, we have the following estimates: 

m s+e °(u)+M Uo+eo {u) + 0{u,uo) + 0(5,u) < C(m ,uo,e ). 

We also have the Sobolev inequalities: 

HWlHSu.u) ^ C(m ,u ,e )(\\f ' <j)\\ L 2 {s ^ u) + \\<p\\ L 2 {s ^ u) ), 

2 

\\<f>\\u°°{S a ,u) < C(m ,u ,e )Y^ \\f l( t>\\L^Su, u )- 

i=0 

For the sake of simplicity, we will suppress the lower index < 3 in the norms, e.g. 
= M<%. In the remaining part of the current section, C refers to various constants 
depending only on mo, u$ and eo and A < B refers to A < CB. 

Since we expect the solution is close to some Schwarzschild space-time, we have to 
derive more precise energy estimates on M £0 . Combined with Rendall's result |14j and 
usual bootstrap arguments, our estimates will be good enough to prove existence. We 
remark that Luk's result [13] is extremely convenient to use at this stage: we can avoid a 
long bootstrap argument and we can take the existence for granted. 

We start to prove Theorem [2] by deriving the estimates only for , and 0k for 
k = 0, 1, 2, 3. In this case, we can take K = 5. The estimates for higher order derivatives 
can be obtained by a routine induction argument. Since it is similar to the higher order 
energy estimates derived in Section [3j we will omit the proof. 

We first prove the following lemma which claims that the connection coefficients can be 
controlled by the initial data and the curvature components. 

Lemma 5. 0(u, u) < sup 5 <„,<„(^(u', u ) + @ u (y/)) + sup Uo < u ,< u (0(6, u') + &*{u')). 
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Proof. The proof is once again by integrating the null structure equations (2.1)-(2.8). At 
the beginning, we remark that for a tensor <f>, we have 

3 3 3 

i=0 i=0 i=0 

The can be derived by the bounds for Q and cHogfi = (ry + rj)/2. Similar inequalities also 
hold for L 4 (S U)U ) and L°° norms. 



We start with the bounds on x- I n view of (2.1), we have 



D\x'\ 2 + 2ntr X \x'\ 2 <2|xlM. 

We then use the L°° bounds for £1 and trx in Proposition[2]as well as the Sobolev inequalities 
to derive 

\\x \\l^(Su,u) ^ llxlU-^) + / \\ a h°°(s u , yU ) 

2 



<E (\\r i nv>is^)+ rw'aWv^S) 

i=0 v J5 7 



We now turn to the derivatives of x- We first rewrite (2.1 ) as Dx' = ^ 2 |x'| 2 ^ ~~ a an d then 
commute with ^7 to derive 



a. 



The last term ^/ a can be bounded by ffl~{u). For the remaining terms, each of them 
contains a factor |x'|- We can bound the other factors by C and bound |x'(m>^)| by 
@(8,u)+M^{u). Thus, we obtain 

\wnLoo {S ^ u) <mu)+^(u). 

The estimates on higher order derivatives are similar. For second and third derivatives of 
X, we have 

fc-l k 

Df k x = £ry (nx) • f^'x' + e^ 2 ) • r k ~Wit - ? k a. 

i=0 i=0 

When k = 2, for the first and the second terms, we can use Holder's inequality by placing 
the factors with highest order derivatives in L 4 (5 njU ) and the others in L°° . In this way, 
for each product, one factor is bounded by 0(8, u) + &-(u) and the other one is bounded 
by C. The last term y k a is also bounded by M-(u). When k = 3, there are two cases: if 
there are factors containing third order derivatives, we can bound it in L 2 (S UjU ) and bound 
the other factors in L°°; if not, we bound the factors containing second order derivatives 
in L 4 (S U! u) and bound the other in L°°. Therefore, we obtain 

iiy'x'iiL^) + \\? 3 x'\\li ( s u _, u ) < mu)+^{u). 
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For x, in a similar manner, we can use (2.3) to derive 



IxlU- + \\fxh~ + Wf xW lhs^ u ) + \\f xh^ Sst , u ) < 0(u,u o ) + M u {u). 



For 77 and 77, we consider (2.5) and (2.6) as a coupled system and we derive 



IMIl-(S„, u ) < \\vh«>(S,, n ) + / IIZZll^^'.J + I \\P\\L<*>(S u , t J, 

Jo Jo 
pu pu 

\\vh°°(Su,u) ^ ll2lU°°(S!b«o) + / Wvh^s^,) + / 

J UQ J UO 

We substitute the second inequality to the first one and take snp s<u / <u on both sides, we 
have 



pu 

sup NUoo (5 u j< \\v\\l°°(S 5 ^ + sup ) + / sup 

Cu'<u o<u'<u Jug 5<u'<u 



S< 



+ / SUP ||£|| L oo (S ) + 
>uo o<u'<u J5 

Thanks to Gronwall and Sobolev inequalities, we have 



pu 

sup \\vh°°(s u , j ^ NIIl-(s 5 ,„) + sup HglU-^^) + / sup 

<5<u'<m <5<«'<u Jun 8<u'<u 



+ 



uo o^.u'^.u 

LOO{Su , j < 0(6, u) + sup 0(u', n ) + sup #"(«') + ^(u) 



/(5 5<u'<M <$<«'<« 

We then define &(u,u) and rewrite the above estimates as 

sup ||r?||z,oc.(„/ )U ) < &(u,u) = sup (0(u',u o ) + ^ a (yf)) + sup (0(6, u') + ^(u')). 

8<u'<u 8<u'<u uq<u'<u 



We turn to the derivatives of 77 and r\. Recall that we can commute derivatives with (2.5) 



and (2.6) to derive (3.14) and (3.16). When k = 1, we have 



pu pu 

Wfvh^s^) Z \\fv\\v»{s 6 ,u) + / Wfvh^s^ + / liy/3IU«(%,j + ^fe,"), 
\\Tv\\L°°( S&n ) < \\?Rh~(s a , U0 ) + / \\Tv\\l~(s &u ,) + / liy^lU-(% u o + 

By a similar argument as above, we obtain 

Similar arguments also apply to the cases when k = 2, 3. We finally obtain 

\\f 2 (v,v)\\LHs^), \\f 3 (v,R)\\Li(s^) < &(u,u). 



For ui and w, the estimates rely on (2.7) and (2.8). We rewrite (12.71) as 



£(w - w mo ) = (fi 2 - J& ) (2(7/, 77) " M 2 " P) + - Pmo)- 



(4.11) 
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In view of the facts that DlogQ = u, |0| + |0 _1 | < C, we have 

\\M ~ ttmo\\L°°(S UlV ) < II lo g^ - l °S^m \\L^(Su,u) 
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< 



\\0L- HUnoU-lS^) £ sup ||w- W m J| LO o (s ). 

UQ UQ<U'<.U 



We can integrate (4.11) to derive 



\u-UL mo \\ L oo( S}ltU ) < 11^-^11^(5^)+ / - n mo \\ LO c {s ^ u) + M-(u) 



< 



'0 uo<u'<u 

Therefore, Gronwall's inequality yields 



sup^ || w -ui mo \\ L oo is ^ ul) +M^(u). 



sup \\U) -WmolU-tS ,) £ sup 

UQ<U'<U U0<u'<U 



As a byproduct, we also have 



|O-O W0 |U~(S s ,„) < sup 



We then rewrite (2.8) as 



D(u -UJ mo ) = (n 2 - 2 lfl )(2(r ? , rf) - \rj\ 2 - p) + fi^(p - p mo ). 
We can integrate this equations and use the bounds for | fi — O mo | to derive 
llw-w^lUoc^) < sup (tfftuO+^V)) 

We turn to the derivatives of to and lo. Since w mo and w mo are constants on Su,u, applying 
J7 k on (2.7) and (2.8) will kill those terms. We also notice that p mo is constant on S u u 



and ddogTT= (t] + rff/2 has already been bounded by We then take derivatives 



k-2 



Df k u> = rnVx) ■ f k - 2 - 1 ^ + f k (n 2 (2( V ,r,) - \rf - p)), 



i=0 
k-2 



Df k u> = Y,Tf{^X) ■ y k - 2 - 1 ^ + f k (n 2 (2(r,,r,) - \rj\ 2 - p)). 



i=0 



For k = 1,2,3, we can integrate as before to derive 

Wifaudhoo^) + liy 2 (^^)ll^(5 M , u ) + \\? 3 (u,l;)\\ l2(s ^ u) < &{u,u). 



Finally, we consider try/ and try/. We rewrite (2.2) as 



1 



D(trx' - try/ ) = -~ (try + try mo )(Otrx' - O mo try/ ) - \x 
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We can bound |Sltrx — ^m trXm | < |fi — f2 mo | + |trx' — trx^J- We then integrate the 
above equation to derive 



||trx / -trx / mo ||L-(5„,„) ^WW -toXmo\\L™{Ss,u)+ / H tr x' ~ tr xL„ IU»(S U /,J 



<5 

+ sup (0(6,u')+& u (u')). 

uo<u'<u 

Thanks to Gronwall's inequality, we have 

||trx , -tr X / mo IU-(S t ., u )^ sup 

In a similar manner, we have the following bound for try/: 

lltrx 7 - trx' IU- (Su , tt ) S n ) + + sup u') + 

«0<u'<u 



We turn to the derivatives of trx' and trx'. We apply to (2.2) to derive 



Ddtrx' = -trx^(^trx') - fi{\x\ ) = -Otrxdtrx' - trytrx'd^ - ${\x\ )■ 

The last two terms can be bounded by &{u,u) and the first term can be absorbed by 
Gronwall's inequality. Thus, we obtain 

\\$to^\\L°°(u,u) % ^>{u,u). 
For higher order derivatives, we have 

Df fc trx' = -Otrxy fc trx' + Lo.t., 

where l.o.t. can be bounded directly by &(u,u). Thanks to Gronwall's inequality, we 
obtain 

\\f Wll L 4 (u,u) L 2 (u,u) ^{Ui U). 

Similarly, we have the following estimates for trx': 

L°° (u,u) + \\f trx'H 

Finally, using the bound of y l (f2 — £l mo ) in terms of &(u,u) for % = 0,1,2,3, the above 
estimates for y 1 (trx' — trx' mo ) and y * (trx' — ^x' mQ ) also hold for y 1 (trx — trx mo ) an d 
y*(trx - trx„J, or f l (Qtrx - ^m tr Xm ) and y l ( fitr X ~ ^m trx mn ), for i = 0, 1,2,3. 

Therefore, we have completed the proof of Lemma [5} I 



We rewrite (2.24) and (2.25) in terms of the renormalized quantities such as p 



Recall that p mo satisfies the following equations 

3 3 

D Pm + 2^mo tr XmoPm = 0, Dpm + « ^ m tr X mQ Pm = 0. 
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We can subtract these two equations from (2.24) and (2.25) to derive 



D (P ~ Pm ) + T^trx - O mo trx mo )p + -ft mo trx mo (p - p mo ) 



3 3 
D{P ~ Pmo) + 2^ ntT K - n "io tr X mo )p + 2°™o tr X m JP - Pmo, 



-m Q 



Because p mo is constant on each S UtU , we can also rewrite (2.22) and (2.23) as 



(4.12) 



(4.13) 



DJ3 + -QtrxP -nx-p + uiP- m(p - Pmo) + *$o + 3VP + 3V + 2* • = 0, (4.14) 
DP + ^ntT X P -Slx-P + u§_ + n{$(p -Pmo)- *fa + 377P - 3* W -2x-P} = 0. (4.15) 



equations (2.18) 



We call (4.12), (4.13), (4.14 



and (4.15) together with the following six original Bianchi 



2.19), (2.20), (2.21), (2.26) and (2.27) the renormalized Bianchi equa 



Hons. In a similar manner as we derived energy inequalities (3.44), (3.45), (3.46) and 



(3.47), we have the following energy inequalities 



E 



\R\ 2 + 



E 



I-RI 



< 



E 



I-RI + 



RE{a,l3,p-p mo ,a,P}' ' w 



E 



I-RI 



E 



lor-, 



fill 

Re{P,p- Pmo ,<7,f3,<x} ~ S 



R£{a,P,p-pm ,<r,P} 
R£{f),P—Pm ,<r,P,a} 
« 1 ,« 2 6{ a ./3.P-Pm ,CT,/3,Qi} 

+|p| |(fitrx - n mo ticxmo){p ~ Pm ) + X ■ " + V ' P + X ■ 9L + V ■ §}■ 
where M UjU refers to the region corresponding to 5 < yf < u, uq < v! < u. We also have 

'*E»|2 



E 



R£{a,f3,p-pm 



7i R\ 2 + 



K ■■■■ ■(/ 



E 



fill 

Re{/3,p-p mo , ( 7,^a}- / ^ 

\f i R\ 2 + J2 I i^i 2 



W R\ 



Re{a,/3,p-p mo ,a,§) - ~«0 
i-1 



R&{/3,p-p mo ,cr,/3,a} 



+ 



E 



Be{/3,p-p mo ,<T,/3,Q:} 



Y,w j f{^x)-y i ^ 3 R-y l R+y j y{^x)-y i ^ j R-y i R\ 



3=0 
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i-1 

+ W jK ■ T'^R ■ f' l R + f j K ■ f^^R ■ f' l R\ 

3=0 

+ \T((nr) • ^) • y^ 2 | + Hiy^mrx - n mo tv Xmo ) • f VI 
+ Hirx ■ y + f \ ■ TP + ft ■ Tql + Tv ■ ?% 

for i = 1,2,3. 

We remark that in terms St\ ■ &2, the term a ■ a do not appear. Thus, we can regard 
&i ■ &2 as either Ri ■ R2 or ij^ • R 2 . We then have the following estimates: 




Therefore, by Cauchy-Schwartz inequality, we have 

3 ~~ 3 ~ u ~ 

£// iTinr-Rj-TRii^ J\TRi\ 2 + \? l R2\ 2 )- 

Similarly, we have 

3 3 

\T(nr. r^-Tr,]^ T [ (ir^il 2 + ir« 2 | 2 ). 

i=0 JJ ^,u i=0 Js J ci, 

We turn to the following terms 

2 

2 — Q ^ Mu^ u 
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They are treated in the same way as above by virtue of the bounds on K. They are 
eventually bounded by 

3 



/ / W l R\ 2 + l W l R? 

Ju Jc% J 8 Jc_l, 



Finally, we consider the following terms: 

E // \p\\f\^ tT X-^mo^Xmo) -y\p-Pmo) 



i=0 uJm vl,^ 



According to Lemma [5j they are bounded by 

snps<u<s +eo ^(u, u Q ) 2 + sup 0(5, u) 2 

uo<u<—l—u 



V / sup / \TR\ 2 + / sup / |V7\R| 2 

u u <u" <u' J C-„ J 8 8<u"<u' Jc u , 



i=0 



We define 



r>^r a oi uq<u'<u • n JC,,i 

f{u) = £ sup £ I 

D^r a „ a ^ S<u'<uZ7. JC„ 



Therefore, the above estimates can be summarized as 
£(u) + T(u) < £(u ) + T(5) 

/ £[u') + I F(v! 

Juq J 8 

Thanks to Gronwall's inequality, we have proved: 
Proposition 3. If we have the following smallness on initial data 

^ 1+5 (u ),M~ 1 ' 5 (5), 0(u, u ), 6(8, u) < e, 
then for (u, u) E [5,8 + eq] x [no, uo + eq], we have 

^(u),St(u),@(u,u),@(u,u) < e. 



The smallness of the data is of course guaranteed by Proposition [TJ 
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fit 



Equipped with Proposition |3j we are able to control C norms. More precisely, we shall 
prove that 

b - 9m \\cs(M eo , grno ) £ e. (4.16) 

We start with the C° norms. We write g in canonical double null coordinates as follows 
g = -2Q 2 (du ® du + du ® du) + $ AB {dB A - b A du) (d6 B - b B du). 

Since P>(0 2 ) = 2ft 2 u;, 1^ = 2tt X and £>& = 4ft 2 C, by virtue of the L°° bounds of w - w mQ , 
X — X mQ an< i C i n Proposition |5|as well as the Sobolev inequalities, we have 



\9 - 9m \\c°(M e 



< 



For C 1 norms of g — g mo , we use V mo (g — g mo ) = (V — V mo )g. We also need the L°° 
bounds on V — V mo , which can be obtained by T — T mo and P — P mQ where Y refers to the 
null connection coefficients and P refers to the Christoffel symbols of fi. We have already 
derived 

|r-r mo | ^ e 

thanks to Proposition [3] and the Sobolev inequalities. To estimate P —f mo , we need the 
following propagation equation 

D(r -T mo ) C AB = \f D (f a^X)bd + ? b (Vx)ad -f d (Vx)ab). 

Since we have already obtained the L°° bounds on ^(0%), by directly integrating the 
above equation, we obtain |P — f mo \ < e. Finally, we have 



9mo\\c^-(M tQ ,g mo ) 



< €. 



For C 2 bounds on g — g mo , we write 

V mo (5 - 9m ) = (V 2 - V 2 mo )g = V((V - V mo )g) + (V - V mo )V mo <?. 

The last term (V — V mo )V mo g has already been controlled. Thus, we need the L°° bounds 
of the following quantities: 

£>(r-r mo ), D(F-f mo ), £(r-r mo ), ^(p-p mo ), f(r-r mo ), f(£-f mo ). 

The estimates for the first four quantities can be obtained using the null structure equations 
provided we have L°° bounds of all first derivatives of null curvature components, but we 
do not have control on 

D(ui-u) mo ), D(ui-ui mo ); 

The estimates for the fifth quantity are obtained directly by Proposition [3j The bounds for 
the last quantity are obtained also by the propagation equation of P — P mQ provided we 
have L°° bound of 

y 2 (%). 
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The strategy is clear now: we have to control [|(-D(w — oj mo )\\L°°-, — ^. mo )\\L°° and 

|y 2 (J7x))l|L 00 - For this purpose, we first claim that 

\\a,P,p- Pm ,cr, f3,a\\ L oo < e. 



This can be obtained by using a variation of (3.4) and (3.5) by setting the constants 
depending on mo, eo, uq and the fact that ||y - ck 1 1 x 4 (s , <5 „) ^5 6 which is stated in Proposition 



1| To apply (|3.4|) and (|3.5|), besides &2,M.2 ^ e > we al so need ||^ Da^iprQs < e and 



D(P,p - Pm o ,0,P,a)hj(c u ) ~ e - For A P ~ Prri , ° and j3 this is done by simply 



taking y to both sides of the renormalized Bianchi equations. For 



||y Da \\r,2(ri ), they are obtained using a variation of (3.48) and (3.49) by replacing Lie 
derivatives by . 



fDa\ 



L 2 (C U ) 



and 



We then claim 



J oo ^ E. 



The proof is similar to the above, we first apply y to (2.22) and (2.21) and then apply 



(3.5) and Sobolev inequalities on S UjU . This yields ||y (/?, /3)||l°° IS e - As a consequence, 
we have 

\\D(uj -aj mo ),D(u - w mo ) 



In fact, we can apply D to (2.8) to derive a propagation equation for D(uj — oo mo ) which 
looks like 

DD(oj - uj ) = D(p - p Q ) + l.o.t. = f/3 + l.o.t.. 
We then integrate this equation to conclude. Similarly, we have estimates for D(co — w TOo ). 

To prove ||y 2 (Slx)ll-L°° ^5 e ) using bounds of \\y s (3\\ L 2( C j given by Proposition 3 we 



couple (2.11) with (2.4) to derive bounds of ||y xIIl 2 (c ) an d \\f ^ t x\\l 2 



Therefore, we have proved 



\9 - 9m \\c 2 (M eo ,g mo ) 



< 



e. 



For higher order estimates on C k norms, as we remarked before, we can simply use a 
routine induction argument and we omit the proof. Therefore, we have completed the proof 
of Theorem [2j I 



5. Gluing Construction 



5.1. Preparation. We summarize some key properties of the metric g constructed in the 
previous sections: 
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(1) In the region (u, u) £ [uq,0] x [no,0], g coincides with the standard Minkowski 
metric; 

(2) In the region where u E [0, 5] and u £ [no, —1 — u\, the metric g is constructed by 



Christodoulou in his work [4J. In addition, we have imposed condition (2.30) with 
|i4q| > 2mo, so that on the incoming cone Cg, the estimates in Section 4.2 hold. 
For a sufficiently small 8, tix(8, uq) > and hence Ss )Uo is not a trapped surface, 
and if 2mo > 1, according to Theorem 17.1 in |4J, Ss~i~s is a trapped surface. 

We remark that, for < rriQ < ^ , it is straightforward from the proof of [1] that 
g actually exists up to u + n = — 2mo and Ss,-2m -6 is a trapped surface, once we 
choose a sufficiently small 8. 
(3) In the region (u, u) £ [5, 8 + eo] x [no, no + eo], the metric g is e-close (in C fc_3 norms) 
to the Schwarzschild metric g mo . 



We fix a sphere Si = Ss, Ul near Ss, Uo by choosing ui close to uq in such a way that 
trx(<5, ui) > and trXmo(<5> ^l) > 0. We emphasize that the choice of Si does not depend 
on 8 if 5 is sufficiently small. Because g mo is the Schwarzschild metric with mass mo, we 
can choose the time function t of g mo in the Boyer-Lindquist coordinates (we regard the 
Schwarzschild metric as a member in the Kerr family) such that t(S,ui) = 0. We also 
regard t as a smooth function in the u-u plane (u, u) £ [8,5 + eo] x [no, no + eo]. 

Now in the region (u,u) £ [5,5 + eo] x [no, no + eo], we choose a three dimensional 
hyper surface 

H = ^ Su,u, 

t(u,u)=0 

and we use S2 to denote the intersection of H with C Uo LiCj +e . We define the following two 
radii r\ = r|5 x and r2 = r|g 2 where r(n, n) = r mo (u, u) is the radius of S M)U in Schwarzschild 
metric 5 mo . By construction, 2mo < ri < r2 and il is space-like with respect to g mo . Thus, 
if e is sufficiently small, H is also space-like with respect to the metric g. 

Let (j mo ,fc mo ) be the induced metric and the second fundamental form of H as a sub- 
manifold in the Schwarzschild space-time with mass rriQ. Therefore, for r\ < r < t2, we 
have 



and 



g mo = (l- 2 -^)-V + rV s2 , 



Let (0, k) be the metric and the second fundamental form of H induced by g, since g is 
e-close to the Schwarzschild metric g mo in C fc_3 norms, we have 

and 

\m\c^-HH,g mo ) ~ e - 
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We recall the form of the Kerr metric in the Boyer-Lindquist coordinates 

2mr 2 2mrasm 2 9 r 2 + a 2 cos 2 9 2 

9m,(0,0,a) =(-1 + o , „ 2 ^2fl) d * ~ .2 ■ „2 ™ c 2 fl dtd< ^ + ,.2 _ w, -L „2 dr 



+ (r 2 + a 2 cos 2 0) 2 d0 2 + sin 2 9(r 2 + a 2 + ^t'^j V- 
v y v r 2 + a 2 cos 2 6* 

where m > a > 0. The lower index (0, 0, a) specifies the angular momentum vector: if we 

use the spherical polar coordinates x 1 = r cost/? sin x 2 = rsin^sin^ and x 3 = rcos6>, 

then the vector (0, 0, a) is the same as the axis (of the rotation of the Kerr black hole) and 

the norm of (0, 0, a) is the angular momentum a. For an arbitrary vector a = (ai, 02, 03) £ 

M 3 with |a| = a and an isometry Q a G £0(3) mapping a to (0,0, a), we can also define a 

family of Kerr metric g m ,a = (idiR x £2 a )* 9 m ,(o,o,a) where ids is the identity map of t axis. 

We remark that this definition does not depend on the choice of ft a . 

We choose the slice H to be t = in the above Kerr space-time and we use (x 1 , x 2 , x 3 ) 
and (r,9,(p) as coordinates on H. Let (g m ,(o,o,a)) &m,(o,o,a)) be the induced metric and the 
second fundamental form of H, thus, we have 



r 2 + a 2 cos 2 6> 2 2 2 2 „n ,„2 / 2 2 2mra 2 sin# , . 

flW»,(0,0,a) = -2 o T^ dr + r + a cos 61 d6 * + ( r + a + 2,2 2fl) Sm 

,y ' ' ' r z — 2mr + or r z + a cos z # 

and 



V(0.M = ^ + fl2+ W^ WaWfl^ " O d ^ 



2m 2 r 2 a 3 sin 26* sin 2 , „ , , 

H r-j- — 2 2 .„ d6>d y)- 

(H + cos^ 6*)^ 

For an arbitrary vector a = (01,02,03) G M 3 with |a| = o, we can also define the cor- 
responding of initial data set by g m ^ = tt* a g m ,(o,o,a) and &m,a = ^ a &m,(o,o,a)- Tne y also 
correspond to the t = slices. 

We now recall the definition for the constraint map that is, 

$(s,tt) = (W(S,7f),diV7f), 

where 7f = k - trkg and W(g,Tr) = R(g) + \{t?Tt) 2 - |7f| 2 . Let D$>*- be the formal L 2 - 
adjoint of the linearized operator of at (g, tt). When (g, tt) = (g mo , 0), we use K to denote 
the null space KerL>$^_ Q y We will explicitly write down K. In (ir^x 2 ,^ 3 ) coordinates, 

we use = 1,2,3) to denote the following vectors fields 

^ 2<5 3 <9 ^ o(9 , (9 _ 1 <9 2 ^ 
S2i = x x - — , tt 2 = x x - — , S2 3 = x x - — , 

OX3 OX2 OXi OX3 OX2 OX\ 

and in (r, 6, ip) coordinates, they can be written as 

Qi = —simpde — costpcotOdp, 2 = cosipde — simp cot 98^,, O3 = d v . (5.1) 
Then we have the following description on K: 
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Lemma 6. K is a 4 dimensional real vector space. In fact, K is spanned by the following 
four elements: ((1 - ^)3, 0),(0, J2i), (0,O 2 ) and (0,fi 3 ). 



Proof. We first recall the form of the linear operator D&(g ^ (see also Section 2.3). For a 
smooth function / and a vector field X, we have -^(f , X) = DH*_ ^(f)+Ddiv^- ^(X), 
where 

Dny } (f) = {{L*f)ij + (tr S 7f • Tfy - 2**0* j)f, (trgif • - 2%)/), 



Ddiv 



U;.7r)P0 = 2^ X7ri 3 + ^kX k lTij - (Xi(VkTT k j + XjVjfc7T* 



- (V m X fc 7r fcm + X k V m Tr mK )gij, -£x9ij), 
and = -(A s f)g + Hess g (/) - /Ric§. In the current case, (g,7f) = (c/ TOo ,0), thus, 



D<5> 



(9m ,0) 



(/,^0 = {L* § f,--Cxg mo ) 



(5.2) 



First of all, we consider the second equation Cx9 



mo 



in (5.2). It amounts to say 



that X is a Killing vector field on constant t slice for g mo . Thus, X is spanned by Ox, 02 
and O3, i.e. the infinitesimal rotations of the Schwarzschild space-time. We remark that 
mo 7^ in this case otherwise we may also have translations of the Minkowski space-time. 



Secondly, we consider the first equation in (5.2): 



L *a f 



-(A- gm J)g mo + Hess 5mo (/) - /Ric 5roo = 0. 



Taking the trace, we obtain Ag mQ / = 0. So the above equation is equivalent to 

Hess s m „ (/) = f Ric S mo ■ 



(5.3) 



We show that / depends only on r. We consider the (d r , S^-component of (5.3), that is, 
drdef—hdef = 0. Thus, def = ar for some aei We then turn to the (d r , (? r )-component 



of (5.3), it can be written as 



m 



37 + ^(1 



2mo , 



2m _ 2m , 



We can take derivative dg on both sides and then substitute dgf = ar to the above equation. 
Those operations yield = — 2 °T' ■ Therefore, a = since mo 7^ 0. We conclude that 
def = 0. We then proceed in the same manner to show d^f = 0. Hence, / depends only 



on r. To obtain the exact form of /, we consider the (dg, (%)-component of (5.3), that is 

m 



d 2 e f + r(l-^)d r f = ^f. 



Since / is simply a function of r, we have ^j- 



[I 



2mp • 



. Therefore, f(r 

It is straightforward to verify that f(r) = a(l 
indeed a solution of (5.3) for all oel. Hence, we complete the proof of the lemma 



ail - 5^)2 for some a £ 



2moU i s 
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5.2. Gluing Construction. We choose a sufficiently large k sufficient in Theorem [T] so 



that we can use the Local Deformation Lemma in Section 2.3 To run the gluing construc- 
tion, we take <ft = <fi{r) a smooth cut-off function on the three slice H such that (ft = 1 near 
Si and (ft = near S2 ■ We define a candidate for the final metric as follows: 

(<?, k) = ((ftg + (1 - (ft)g m , a , (ftk + (1 - <ft)k m ,s), 



At this stage, (g, h) may not sat isfy the constraint equations (1.1) and we will use the 



techniques introduced in Section 2.3 to deform (g, k) to be a solution of the constraint 
equations. 

We choose two free parameters m 6 IR>o and a£R 3 , so that 

\m — mo\ + |a| < Coe, 

where Co is large positive constant to be determined later. The constant Co is independent 
of the choice of e, m and a. 

By the construction, it is easy to see that 



\g - 9m \\ c k~4, a < e, ||fc|| c fc-5, Q < e, 



for some < a < 1 and p is a weight function decays near dH = S\ U £2 (see Section 2.3 
for definitions). We remark that the above norms are computed with respect to g mo . Since 
both (g, k) and {g m ,a-, k m ,a) solve the constraint equations, it is easy to see that 



k — 6,a 



c _r 



< 



where tt = g — trgkg. 



We now follow the procedure in pQ and [2] to deform (g, n). Let C be a bump function 
on H which is compactly supported between Si and S2, we define (K to be the following 
set vector space 

(K = {((fXX)\(f,x)eK}. 

It is the obstruction space for the gluing for data close to the Schwarzschild data. We 
remark that if £ is sufficiently close to constant 1 and if the data (g, tt) is sufficiently to 
those of the Schwarzschild data (g mo ,0), then the restriction of the operator DQ*, , on 
the orthogonal complement of (,K is still injective. In this case, we can still apply the 



Local Deformation Theorem (Theorem 2 of [2] or Section 2.3) to deform the image into QK 
instead of K. We also notice that, in this case, if we use V = L 2 (C(H)) x L 2 ((S(H)) to 
denote the pair of square integrable symmetric two tensors on H, thus, V = K © ((K)-*-. 
After applying the Local Deformation Lemma, there will be no component in ((K) 1 - . We 
also fix a weight function p which behaves as p ~ d(-,dH) N for sufficiently large N near the 
boundary. According to the Local Deformation Theorem , if e is sufficiently small, there 
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exists a pair (h,uj) G S k ~^ a (H) x S k ~ 5 ' a (H) such that <f>(g + h,n + ui) G = 

near dH = S\ U S^- Moreover, (h,uj) satisfies the following estimates 

||(/i,w)|| c , fe -5, Q < e. 
p 1 

We also remark that, instead of the polynomial weight p defined above, we can also use 
an exponential weight p, namely p ~ e - ^ near the boundary, to solve (h,uj) in the smooth 
class instead of Holder classes (because 3>(g,7r) is strictly supported in the interior of H). 
One can refer to pQ and [2] for more details. 

We are going to argue that we can choose suitable parameters (m, a) so that <&(g + h, 7r + 
ui) = 0. According to the above discussion, in order to show this, it suffices to show the 
projection of &(g + h,n + ui) to K vanishes, i.e. the following projection map X hits zero: 

X : B Coe -»• K 4 , (m,a) i-> (X 0l X 1 ,X 2 ,X 3 ), 

where -Bcoe = {( m ) a ) £ M 4 ||m — mo| + |a| < Coe}, 



and 



for % = 1, 2, 3. 



2o(m,a)= / (1 - — )^H(g + h,n + ui)dp 
Jh r 

Xi(m, a) = / divg +/l (fp + w)jJ^d// 5 . 

JH 



We first analyze Xq. We consider the Taylor expansion of % in Banach space C k 5,a 
near the point (g, fr), thus, we have 



Z (m,a) = / (l-^)i(^(g, 7 r) + Z)^ ( ^ ) (/ i ,^) + 0(||(/ i , w )||^_ 5 , a ))d^ mo . 



i - i - 1 n \ .</. /i ) -r u n, (/ ._, i //. ^ ; -r || \ " ■ -~ 

I H 

We observe that the second and the last terms in the above integral are of size 0(e 2 ) with 
respect to our small parameter e. This is obvious for the last term. For the second term, 
recall that (1 — ^r 1 )^ G KeiDHI- n , and (h, uj) vanishes near dH, therefore, 

iH ^- 2 ^Dn iSmofl) (h,u J )dp- 9mo =o. 

We then subtract the above quantity from the second term to retrieve one more e as follows 

^- 2 ^) 1 'Dn (§ ^(h,u)dp- gmo = [ (i-^)h(Dn i ^ ) -Dn i§m(v0) )(h 1 u J )dp- 9mo 



< 



9m ,Tr)\\c4( h ,u)\\c2 = 0(e 



2\ 



for k > 7. Finally, we consider the first term in Xq. Because |a| < Coe, we will make use of 
the following key observation: 

5 mia= (l-^)-Mr 2 + r 2 d4 2 +0(6 2 ). 
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This is obvious according to the formula for g m ^ in last section. We also write 

g = (1 _^)-V + r 2 dff 2 2+/lt 

where h e is a two tensor and \h e \ < e. Therefore, we have 

f (i_^0)J Wj#)d/i =f R(f(r)dr 2 + r 2 da 2 s2 )dn- go 
Jh r Jh 

+ Jjl- ^) § ^(^ ,0)(#e,0)d^ mo +0(6 2 ) 

where /(r) = 0(1 — ^r 1 ) -1 + (1 — 4>)(1 — 2 y L )~ 1 - To compute the first term on the right hand 
side, we have to compute the scalar curvature R(f(r)dr 2 + r 2 da 2 ). By direct computation, 

R{f{r)dr 2 + r 2 da 2 ) = 2 r - 2 (rf- 2 (r)f(r) - f-\r) + 1) = -2r" 2 a r (r(/- 1 (r) - 1)). 

By virtue of this formula, we can easily derive 



R(f{r)dr 2 + r 2 da 2 s2 )dfi- go = -8^(r(/~ 1 (0 " = 167r ( m " m o). 

For the second term, we denote 



g A _ f (1 _ ^)h D H M (cj>h e ,0)d» gmo . 

J H 



Recall that (1 — ^r 1 ) 2 e KevDH*^ y <fi = 1 near Si and = near S2, therefore eo can 

be written as an integral on S±, which only depends on h e on S±, and |eo| < e. Finally, we 
summarize the above calculations for Iq as 

X (m, a) = 16vr(m - m ) - e + 0(e 2 ). (5.4) 

For each i = 1, 2, 3, we analyze Zj as follows: 

li(m,a)= / div^^Tr + w^d^ + / (div gm - div gmo )(* + wjj^d^ 

= / div <? mo (* + " W'd^ + 0(e 2 ) 

Recall that div| (J7j) = 0, (g, tt) coincides with (g, fc) on Si and (g, n) coincides with 
(p m) a, 7f m a ) on S*2. We can apply Stokes Theorem to derive 



(m,a)= / (km,ediM{dr) l dn r2 - f kM(d r ) l dfi ri +0(e 2 ), 

JS 9 . J Si 



where d/x r is the volume form on the round sphere with radius r. Let ej = J" s fy?^ (^r)'d/x n 
and it is a constant coming from the integral on the inner sphere S±. Since g is close to g mo , 



60 JUNBIN LI AND PIN YU 

we know that |ej| < e. We emphasize that ik is independent of (m, a) and the constant Co- 
Thus, 

2i(m, a) = / (fe m ,a)o-^i (dr) l dfj, r2 -£i + 0(e 2 ). 
To understand the local behavior of Xj, we study the differential of the map 

J:a^(/ (fc m , a )^i(dr)'d/v 2 , / (fc m , a )^(<9 r yd^ r2 , / (A- m , a )«j^3(9 r )'d/i r2 ), 

%/ S2 ^ 52 ^ S2 

at a = and the parameter m is fixed. 
We first compute 

d f 

dJ | a =o (0,0,1) = — \ t=0 ( (k m){ o,o,t))ij^i{d r ) l dfi r2 ) i=1: 2,3- 
When i = 1, we have 

(fcm,(0,0,t 



T2 



= - d6 cosipcotOk rn ^ fl ^(d v ,d r )\ r =r 2 ■ r^sinM^ = 0. 
Jo Jo 

since fc m ,(o,o,t)(^>) &r) does not depend on 99 and J Q 27r cos </?d</? = 0. When i = 2, similarly, 
we have 

/ {km,(0,0,t))lj^U d r) l dfX r2 =0. 
J S2 

When i = 3, in view of the exact formula for fc m ,(o,o,a) ™ previous section, the drd(p 
component of k m ,(o,o,t) is of the form + 0(t 2 ) for sufficiently small t. Thus, 

r f'TT /»27T 

/ (k m ,(o,o,t)h n U d r) l dVr 2 = d6 2mtsm9dip + 0{t 2 ) 
Js 2 Jo Jo 

= 8-irmt + 0(t 2 ). 

Therefore, we finally have 

d 

\t=o { 

>s 2 

or equivalently dJ | a= o (0, 0, 1) = (0, 0, Sum). 



=0 ( / {km,(0,0,t))lj^Udr) l d^ r2 ) = ^TTTTl, 
J So 



I (T 

To compute dJ | a= o (0, 1, 0), we use the SO(3) symmetry on a. We take R = I 1 

\0 -1 0, 

and use dR to denote its differential. Apparently, we have 

dR : (fii, n 2 , fi 3 ) ^ -«3, ^2). 
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Therefore, we have 

d f 

dJ | a =0 (0,1,0) = — | t=0 ( (R*k m){0;0 ,t))lM(dr) l dlJ,r 2 )i=l,2,3 



at 

( / (km,(0,0,t))lj 



dt 



S-2 



= (0, -8vrm,0). 
l\ 

Similarly, we take R = | 10 and we gave 

-1 0/ 

dR ■. (nx,tt 2 ,n 3 ) ^ (-ns,n 2 M. 

As a result, we obtain 

dJ | a=0 (0, 1,0) = (-8vrm, 0, 0). 

Combining all those computations on dJ | a =o and the fact that J(0) = 0, for (m, a) G Bc ot , 
we have 

/-l 0\ 
J(a) = 8vrm -1 a + 0(e 2 ) 

V° 01/ 

f-l 0\ 

= 87rm -10 a + 0(e 2 ). 

V °V 

We now understand the local behavior of X near (mo, 0): 

Z(m, a) = (167r(m - m), 87rm (-ai, -a 2 , a 3 )) - (e , e») + 0(e 2 ), (5.5) 
where |ej| < e. 

We now carry out a degree argument as in pQ and [2] to show that Z(m, a) = for some 
(m, a) G Bc oe provided Co is sufficiently large. We define 

Xi(m,a) = (167r(m- mo), 87rmo(-ai, -02,03)) - (eo,£fc), 
and we then choose Co such that 167rCo > 2|eo| and 87rmoCo > 2|efc|. Therefore, we have 

Zi(m + -^ ) -,— (-ei,-e 2 ,e 3 )) = 0, 

107T OTTmo 

and (m + j^, §^(-ei, -e 2 , e 3 )) G -B Coe . Therefore, Xi is a homomorphism from 5 Coe 
to another box containing in IR 4 centered at (eo, €&). We then define 

I(t, m, a) = (167r(m - m ), 87rm (-ai, -a 2 , a 3 )) - (eo, e k ) + tO(e 2 ) 

on [0, 1] x Bce to be an homotopy between Xi and I. For sufficiently small e, ^ X([0, 1] x 
dBce )-, so the degree of X at the value is equal to the degree of X\ at the value 0, which 
is 1. This implies that X(m, a) = for some (m, a) G Bc oe . 
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Therefore, by a suitable choice of (m, a) 6 -Bc e> we finally deform (5, k) to be (55, fcs) 
where gs = g + h and k$ = ^ + ui + §trg g (fr + oj)gs, so that (55, feg) satisfies the constraint 
equations on H, (g s ,k s ) = (g,k) near Si and (gs,ks)) = (s m ,a,L,a) near 5 2 . This 
completes the gluing construction. 



5.3. Proof of the Main Theorem. Recall the choice of H, it is only defined for u > 5 
and it has an inner boundary at H n Cg. We would like to extend H to the interior to 
complete it as a three manifold without boundary. This can be done by choosing a smooth 
function u = f(u), such that 

'/'(¥) = -1, forn<0, 
f'(u) < 0, for < u < 5, 
t(uj(u)) = 0, foru>5. 
In addition, / can be chosen such that trx(u, f(u)) > for all < u < 5. The space-like 



piece in the picture in Section \1.2\ can be viewed as the graph of / in u-u plane. Suppose 
that the curve u = f(u) intersects the "central" line u = u at (u ccn , f(u ccn ) = u C en)- Let 
=f(u),u cen <u<S Su,u which is smooth space-like hypersurface of M with (g, k) as the 
induced metric and the second fundamental form. It is clear from the construction that 
(g,k) = (<%,0) for u cen < u < where 5 is the Euclidean metric on the ball of radius ro 
which is the radius of So,/(o) computed in g. It is also clear that < ro < r±. 

We divide S = M 3 with Cartesian coordinate (x 1 ,^ 2 ,^ 3 ) into four concentric regions 
Em, £0, and as 

S M = {x\\x\ < r }, S c = {x\r < \x\ < n}, 

Y, s = {x\ri < \x\ < r 2 }, T, K = {x\\x\ > r 2 }. 

By construction we have r\ — ro = 0(5) and r 2 — r\ = 0(cq). Thus, Hq is diffeomorphic to 
Sm U Sc and H is diffeomorphic to £5. These diffeomorphisms are realized through the 
Cartesian coordinates in an obvious way. We now define ((/£, ks) on S as follows 

(g,k), on S M US C , 
fe,^s) = < (gs,ks), on S s , 

fcm,a)) on T, K . 

By construction, (g-%, k-%) is smooth on E. It is Minknowski inside, namely, (<?e, = (S, 0) 
on Em; similarly, it is Kerr on the outer region with \m — mo\ + |a| < e. Moreover, 
according to Christodoulou'w work [4], there will be a trapped surface in the future domain 
of dependence of Em U Sc. 

To complete the proof of the Main Theorem, we also have to show the non-existence of 
trapped surfaces on E. In fact, for any two sphere S embedding in E. Let S rs = {x\\x\ = 
r$} be the sphere which is the innermost one of the spheres in the form S r = {x\\x\ = r} 
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containing S. Then S and S rs tangent at some point p. The outer null expansion of S, 
9s = trsfcs + Hs where tr^fes is the trace of k^, computed on S and Hs is the mean 
curvature of S in S with respect to the normal pointing out. By construction 6$ r > 
for all r. Then at the point p, tr sksip) = trs r5 kz(p) since S and S rs tangent at p, and 
Hsip) > Hs r (p) by the maximum principle. So 

Osip) = tr s ~k s (p) + Hsip) > tr Srs ~kx(p) + H Srs (p) = #s rs (p) > 

and then we conclude that S is not trapped at the point p. Therefore, we complete the 
proof of the Main Theorem. 
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